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ABSTRACT 


Two  metal  combustor  system  concepts  employing  the  lithium- 
sulfur  hexafluoride  reactant  combination  are  currently  under 
development-  One  concept  is  the  batch  or  pot-type  combustor  (where 
fuel  and  combustion  products  remain  within  the  combustion  chamber), 
and  the  other  is  the  steady  combustor  (where  fuel  continuously  enters 
and  products  continuously  leave  the  combustion  chamber).  Both 
concepts  are  desirable  for  application  as  a  thermal  energy  source  for 
underwater  power  systems  utilizing  closed  thermodynamic  cycles  since 
the  lithium-sulfur  hexafluoride  reactant  combination  has  a  high  energy 
density  and  no  requirement  for  combustion  product  exhaust  to  the 
ambient  environment.  The  steady  combustor  concept  is  advantageous 
for  long  duration  operation,  where  the  volume  of  the  fuel  inventory 
is  large  in  comparison  to  the  volume  required  to  provide  sufficient 
heat  exchanger  area  and  reaction  volume.  In  this  situation,  the  use 
of  a  pot-type  combustor  requires  the  presence  of  an  excessive 
quantity  of  fuel  in  the  active  zone  of  the  combustor. 

This  report  documents  the  results  of  a  developmental  program 
concerning  steady  metal  combustors-  Two  combustors  are  described  which 
have  demonstrated  operation  for  periods  of  up  to  ten  hours  at  power 
levels  up  to  25  kW.  The  thermodynamic  characteristics  of  the  ternary 
immiscible  liquid  mixture  present  in  the  reaction  chamber  were 
examined  analytically  by  developing  a  theoretical  model  for  the  mixture, 
and  experimentally  by  generating  density  and  solubility  data.  The 
model  permits  prediction  of  system  characteristics  at  operating 
conditions  not  examined  experimentally,  and  is  useful  in  the  analyses 


CHAPTER  I 


INTRODUCTION 

1.1  General  Statement  of  the  Problem 

The  increased  use  of  closed  thermodynamic  heat-power  cycles  for 
mechanical  power  in  underwater  applications  has  motivated  development 
of  thermal  energy  sources.  A  thermal  energy  source  is  a  source  of 
heat  at  an  elevated  temperature ,  and  can  be  provided  by  exothermic 
nuclear  or  chemical  reactions.  A  thermal  energy  source  can  be  combined 
with  various  closed  thermodynamic  heat-power  cycles  (e.g. ,  Rankine, 

Bray ton,  Stirling,  etc.)  to  provide  a  mechanical  or  electrical  power 
system.  The  objective  of  the  present  study  was  to  investigate  both 
analytically  and  experimentally  the  characteristics  of  a  steadily 
operating  metal  combustor  thermal  energy  source. 

Features  which  make  a  thermal  energy  source  attractive  for 
underwater  applications  include  the  following:  a)  the  reactant 
combination  should  have  a  high  energy  density  on  both  a  mass  and 
volumetric  basis;  b)  the  reaction  products  should  be  nongaseous  with 
a  volume  approximately  equal  to  that  of  the  reactants;  c)  the 
reaction  process  should  be  controllable  within  temperature  limits  of 
available  construction  materials;  d)  the  reactants  and  products  should 
be  nontoxic;  e)  the  reactants  and  products  should  be  compatible  with 
the  construction  materials;  f)  the  reactants  should  be  stable  and 
storable  for  extended  periods  of  time;  and  g)  the  reactants  should  be 
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inexpensive.  Features  (a)  and  (b)  are  particularly  important  for  deep 
depth  operation,  where  large  vehicle  volumes  and  the  design  of 
product  exhaust  systems  for  gaseous  combustion  products  become 
problematical. 

Thermal  energy  sources  having  these  features  are  provided  by  the 
general  class  of  chemical  reactions  between  alkali  metals  and  the 
halogen  or  halogenated  gases.  The  particular  member  of  this  class 
considered  in  the  present  study  is  the  reaction  between  lithium  metal 
and  sulfur  hexafluoride  gas,  proposed  by  Pauliukonis  [1],  The  lithium- 
sulfur  hexafluoride  system  was  chosen  due  to  its  extremely  good  mass 
and  volumetric  energy  density  and  the  nontoxic  nature  of  the  oxidant. 

Van  der  Sluijs  [2]  reviewed  the  many  possible  reactant  combinations 
and  concluded  that  the  combination  chosen  for  study  is  the  best  for 
application  as  an  underwater  thermal  energy  source. 

1.2  Description  of  the  Process 

1.2.1  Reaction  Characteristics.  The  reaction  is  conducted  in  a 
combustion  chamber  containing  the  molten  fuel.  The  oxidizer  enters 
the  chamber  as  a  gas  through  an  injector.  The  overall  stoichiometry 
of  the  lithium-sulfur  hexafluoride  reaction  is 

8  Li  (1)  +  SFfi  (g)  -  6  LiF  (1)  +  Li2S  (1) 

The  products  of  the  reaction  are  liquids  for  the  usual  range  of 
operating  conditions,  i.e.,  bath  temperatures  greater  than  approximately 
1075°K.  The  volume  of  the  products  is  approximately  96%  of  the  volume 
of  the  original  fuel  at  a  temperature  of  1200°K. 

Characteristic  of  the  phase  equilibria  between  liquid  alkali 
metals  and  their  salts  is  the  formation  of  an  immiscible  liquid  region. 
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At  normal  operating  temperatures  of  the  combustor,  1090°K-1260°K,  the 
combustor  bath  is  within  this  immiscible  liquid  region,  and  the  bath 
splits  into  two  phases.  The  fuel-rich  phase  has  a  molar  concentration 
of  fuel  greater  than  95%,  and  the  product-rich  phase  has  a  molar 
concentration  of  products  greater  than  99%.  The  density  of  the 
product-rich  phase  is  more  than  three  times  that  of  the  fuel-rich 
phase  and  the  products  settle  under  the  action  of  gravity  to  the  bottom 
of  the  combustion  chamber  where  they  can  be  allowed  to  accumulate 
(pot-type  combustor  concept),  or  can  be  continuously  removed  while 
fuel  is  continuously  added  (steady  combustor  concept).  The  chemical 
energy  released  by  the  reaction  is  transferred  to  the  power  cycle  by 
means  of  heat  exchanger  surfaces  located  at  the  walls  of  th'j 
combustion  chamber. 

1.2.2  Pot-Type  Combustor  Concept.  The  pot-type  combustor  is 
illustrated  schematically  in  Figure  1.  The  combustion  chamber  is 
sized  to  carry  the  entire  fuel  inventory  for  the  particular 
application.  Operation  is  terminated  when  all  the  fuel  is  consumed 
and  the  combustor  is  filled  with  reaction  products.  The  oxidizer  is 
introduced  to  the  combustion  chamber  through  an  injector.  Both  single¬ 
passage  injectors  which  can  be  cooled  to  decrease  the  susceptibility 
of  the  injector  to  corrosive  attack  by  the  oxidizer;  and  coaxial 
injectors,  which  utilize  a  secondary  flow  of  argon  gas  in  a  shroud 
tube  around  the  oxidizer  core  tube  to  prevent  contact  of  the  oxidizer 
with  the  outer  injector  parts,  have  been  reported  in  the  literature. 

In  the  case  where  a  single-passage  injector  is  used,  the  combustion 
chamber  can  be  sealed  and  the  chamber  pressure  will  be  equal  to  the 
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vapor  pressure  of  the  liquid  mixture.  In  the  case  where  a  coaxial 
injector  is  used,  the  combustion  chamber  must  be  vented  to  remove  the 
shroud  argon  flow  and  the  combustion  chamber  pressure  is  controlled 
as  desired.  Thermal  energy  is  transferred  through  the  walls  of  the 
combustion  chamber  to  the  working  fluid  of  the  power  cycle  througn  the 
use  of  a  heat  exchanger  which  is  usually  integral  with  the  combustion 
chamber  walls. 

Several  pot-type  combustor  systems  are  discussed  in  the 
literature.  The  earliest  documented  design  is  described  in  a  United 
States  Patent  by  Percival  [3],  who  proposed  the  use  of  Freon,*  which 
contains  fluorine,  as  the  oxidizer.  Distinguishing  features  of  the 
design  include  submerged  injection  through  a  dispersion  pipe 
containing  multiple  nozzles,  and  energy  removal  from  the  bath  by  a 
direct  contact  heat  exchanger,  where  the  working  fluid  for  the  heat 
engine  cycle  passes  directly  through  the  combustion  chamber. 

Uhlemann,  Spigt  and  Hermans  [A]  have  described  the  development 
of  a  pot-type  metal  combustor  system  employing  lithium  and  sulfur 
hexafluoride  for  use  with  a  Stirling  engine.  Bierman  [51  has  also 
discussed  a  lithium-sulfur  hexafluoride  combustor  system  to  be  used  in 
conjunction  with  a  Stirling  engine.  In  this  system  the  inside  walls 
of  the  reaction  chamber  were  covered  with  a  capillary  wick  which  was 
completely  filled  with  lithium.  The  oxidizer  was  injected  into  the 
gas  phase  region  of  the  chamber  and  reacted  at  the  surface  of  the 
wick.  In  this  manner,  the  reaction  chamber  walls  were  protected  from 
corrosive  attack  by  the  oxidizer.  The  reaction  products  formed 
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droplets  on  the  surface  of  the  wick  and  fell  to  the  bottom  of  the 
reaction  chamber.  The  energy  of  reaction  was  removed  by  a  heat 
exchanger  located  behind  the  wick, 

Mattavi,  Heffner  and  Miklos  [6]  also  discussed  the  use  of  metal 
combustion  systems  in  combination  with  a  Stirling  heat  engine.  Several 
oxidizers  were  investigated  including  sulfur  hexafluoride.  The 
oxidizer  was  injected  above  the  molten  lithium  surface  for  some 
tests;  on  other  occasions  a  submerged  injector  was  used.  The  authors 
presented  system  performance  information  in  terms  of  system  weight 
versus  run  length  and  average  power  level. 

Van  der  Sluijs  [2]  presented  a  discussion  of  a  lithium-sulfur 
hexafluoride  combustion  system  to  be  used  as  the  thermal  energy  source 
for  a  Stirling  engine.  The  lithium  fuel  was  contained  in  a  reactor 
vessel  along  with  a  quantity  of  sodium,  and  the  oxidizer  was  injected 
through  a  submerged  injector.  The  injector  design  consisted  of  a 
orifice  opening  which  could  be  closed  by  a  needle.  The  needle  position 
was  controlled  by  the  oxidizer  supply  pressure  acting  on  a  sealed 
piston,  and  on-off  operation  of  the  injector  could  be  controlled  by 
the  position  of  an  oxidizer  control  valve.  The  heat  transport  to  the 
Stirling  engine  was  accomplished  by  evaporation  of  the  sodium 
contained  in  the  reactor  bath  and  subsequent  condensation  on  the 
engine  heater  tubes.  The  condensed  sodium  was  then  transported  back 
to  the  reactor  bath  using  an  electromagnetic  pump. 

Faeth  and  co-workers  [7,  8]  investigated  pot-type  combustors 
employing  the  lithium-sulfur  hexafluoride  reactant  combination. 
Particular  emphasis  was  placed  on  long  term  operation  (order  of  hours) 
with  highly  variable  thermal  load  requirements.  Both  vented  and 
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unvented  combustor  configurations  were  studied.  Various  injector 
configurations  were  examined  including  cooled  and  uncooled  single¬ 
passage  injectors  and  coaxial  injectors.  Thermal  energy  was  removed 
both  by  heating  air  as  a  working  fluid  and  by  direct  radiation  from 
the  combustion  chamber  walls  to  the  surroundings. 

1.2.3  Steady  Combustor  Concept.  For  long  duration  operation, 
the  volume  of  the  fuel  inventory  can  become  large  in  comparison  to  the 
volume  required  to  provide  sufficient  heat  exchanger  area  and  reaction 
volume.  In  this  situation,  the  use  of  a  pot-type  combustor  requires 
the  presence  of  an  excessive  quantity  of  fuel  in  the  active  zone  of 
the  combustor.  A  large  combustion  chamber  volume  causes  several 
problems:  poor  thermal  response;  excessive  parasitic  heat  loss; 
structural  difficulties  of  large-size  high-temperature  components; 
system  size  and  weight  penalties;  potentially  poor  bath  heat  transfer 
characteristics;  and  potentially  poor  fuel  utilization. 

An  improved  approach  over  the  pot-type  combustor  system  for 
long-duration  operation  involves  the  steady  combustor  concept, 
illustrated  schematically  in  Figure  2.  In  this  case,  the  fuel  and 
oxidizer  enter  the  combustor  and  the  products  are  drawn  off  in  a 
continuous  fashion.  The  size  of  the  combustion  chamber  can  be 
optimized  with  respect  to  heat  exchanger  area  or  reaction  volume 
requirements  —  whichever  is  limiting. 

The  operation  of  the  steady  combustor  involves  the  use  of  the 
immiscibility  characteristics  of  the  fuel  and  product  [9].  In  the 
absence  of  strong  stirring,  the  heavy  product  liquid  settles  to  the 
bottom  of  the  combustor  where  it  is  removed  through  a  trap.  Fuel 
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can  be  added  steadily,  or  intermit tantly  since  a  small  change  in 
product  level  in  the  trap  can  accommodate  relatively  large  changes  in 
combustor  fuel  level  (since  the  product  liquid  density  is  approximately 
three  times  the  density  of  the  fuel  liquid).  Both  single-passage  and 
coaxial  injectors  can  be  used  with  the  system.  The  thermal  energy 
is  removed  by  a  heat  exchanger  installed  at  the  combustion  chamber 
walls.  Details  of  the  particular  steady  combustor  to  be  used  in  the 
course  of  the  present  study  will  be  described  in  Chapter  IV. 

A  search  of  the  literature  revealed  no  previous  studies 
concerning  a  steady  metal  combustor  system  other  than  the  concept 
presented  by  Feath  [9],  Preliminary  work  for  this  study  is  reported 
in  References  [7]  and  [8]. 

1.3  Previous  Related  Studies 

1.3.1  Thermodynamics .  The  performance  of  the  lithium-sulfur 
hexafluoride  system  in  terms  of  energy  release  rate  has  been  discussed 
in  the  literature.  Van  der  Sluijs  [2]  presented  an  estimate  of  the 
thermal  energy  released  by  the  reaction  and  investigated  system 
performance  for  a  pot-type  combustor  system.  Biermann  [5]  presented 
performance  information  for  a  particular  pot-type  combustion  chamber 
in  terms  of  energy  release  rate  per  volume  and  per  surface  area  of  the 
combustion  chamber.  Mattavi,  et  al. ,  [6]  investigated  the  overall 
performance  of  a  pot-type  combustor  system  with  respect  to  total 
system  mass  and  system  volume  as  a  function  of  run  length  and  average 
power  requirements. 

The  results  of  the  performance  studies  discussed  above  are 
generally  limited  in  applicability  to  a  particular  combustor  system. 
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and  are  not  easily  extended  to  the  analyses  of  other  systems.  Another 
limitation  is  the  lack  of  consideration  of  the  thermodynamic 
properties  of  the  liquid  and  gas  phases  present  in  the  combustor. 

The  thermodynamic  information  is  important  in  analyses  of  all  combustor 
systems.  For  example,  knowledge  of  the  fuel  solubility  in  the  product- 
rich  phase  is  required  since  the  fuel  dissolved  in  the  product 
potentially  represents  a  fuel  loss  in  the  steady  combustor  concept 
and  fuel  unavailable  for  combustion  in  the  pot-type  combustor  concept. 

Information  is  available  concerning  the  thermodynamic  properties 
of  a  majority  of  the  reactants  and  products  in  the  lithium-sulfur 
hexafluoride  reaction.  The  thermodynamic  properties  of  sulfur 
hexafluoride,  lithium  fluoride  and  lithium  are  tabulated  in  the  JANAF 
Thermochemical  Tables  [10].  The  densities  of  lithium  fluoride  and 
lithium  as  a  function  of  temperature  are  given  by  Janz  [11]  and  Davison 
[12],  respectively. 

The  information  on  lithium  sulfide  is  less  abundant.  The  heat 
of  formation  of  solid  lithium  sulfide  was  calculated  by  Morris  [13], 
and  determined  experimentally  by  Juza  and  Uphoff  [14],  The  melting 
and  boiling  points  were  determined  experimentally  by  Cunningham,  et  al., 
[15]  and  Mott  [16] .respectively .  The  density  of  solid  lithium 
sulfide  at  a  temperature  of  298°K  has  been  reported  by  Mourlot  [17], 
Zintl,  et  al.,  [18],  and  Beck  [19]. 

Dworkin,  Bronstein,  and  Bredig  [20]  presented  solubility  data 
for  the  immiscible  liquid  region  of  the  lithium-lithium  fluoride 
system.  The  lithium-litnium  sulfide  and  lithium  fluoride-lithium 
sulfide  systems  have  received  less  attention.  Faeth  [21]  investigated 
the  freezing  point  depression  of  pure  lithium  upon  the  addition  of 


lithium  sulfide.  These  findings  were  supplemented  by  the  investigations 
of  Cunningham,  et  al. ,  [15].  Faeth  [21]  also  determined  the  liquidus 
line  representing  the  freezing  point  depression  of  pure  lithium 
fluoride  upon  the  addition  of  lithium  sulfide. 

Various  thermodynamic  models  exist  for  multicomponent  solutions 
and  will  be  discussed  in  Chapter  II.  Avery  [22]  employed  the  van  Laar 
regular  solution  model  for  lithium- lithium  fluoride  and  sodium-sodium 
chloride  binary  solutions.  Nakanishi  [23]  employed  the  Margules 
model  to  compute  the  eutectic  point  of  various  ternary  molten  salt 
mixtures. 


1.3.2  Liquid  Metal  Combustion  Studies.  Several  studies  are 
reported  which  are  concerned  with  the  combustion  of  liquid  metals 
with  halogen  or  halogenated  gases.  Mattavi  [24]  investigated 
analytically  and  experimentally  the  reaction  characteristics  of  a 
chlorine  bubble  moving  through  liquid  sodium.  Little  [25]  examined 
the  reaction  characteristics  of  a  static  molten  lithium  metal  surface 
during  the  first  few  seconds  after  exposure  to  nitrogen,  oxygen, 
dichlorodifluoromethane,  octof luorocyclobutane,  and  sulfur  hexafluoride. 
Avery  [22,  26]  developed  a  correlation  for  the  combustion  length  of  an 
oxidizer  jet  discharging  from  a  choked  submerged  injector  by  using 
data  on  a  chlorine  jet  discharging  into  a  bath  of  sodium  and  sodium 
chloride.  The  results  of  this  study  are  of  interest  in  combustion 
chamber  design,  since  the  oxidizer  jet  must  remain  within  the 
confines  of  the  bath  if  corrosive  damage  to  the  combustion  chamber  is 


to  be  avoided. 
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1.4  Specific  Statement  of  the  Problem 

Previous  studies  of  liquid  metal  combustion  have  been  concerned 
with  general  hardware  development  and  specific  details  of  the 
combustion  process.  No  attempt  has  been  made  to  develop  an  adequate 
thermodynamic  model  for  the  liquid  and  gas  phases  present  in  the 
combustion  chamber.  Essential  design  information,  such  as  combustor 
bath  composition  and  density,  the  cover  gas  composition,  and  system 
performance  are  also  unknown  at  this  time. 

The  present  investigation  considers  the  development  of  a 
thermodynamic  model  for  the  ternary  immiscible  liquid  mixture 
present  in  the  reaction  chamber  of  a  combustor  using  the  lithium- 
sulfur  hexafluoride  reactant  combination.  The  data  required  to  test 
the  model  have  been  generated  using  a  steady  vented  combustor,  which 
provides  steady  state  conditions  for  the  measurements.  Steady  state 
conditions  allow  accurate  performance  measurements  to  be  made,  and 
meaningful  bath  samples  to  be  taken.  The  present  investigation  also 
demonstrates  long  term  operation  of  a  steady  metal  combustion  system, 
which  had  not  been  achieved  previously. 

The  thermodynamic  model  developed  in  the  course  of  this  study 
will  be  useful  for  the  analyses  of  other  combustion  processes  and 
combustor  concepts  involving  the  lithium-sulfur  hexafluoride  reactant 
combination.  An  example  of  such  an  application  is  the  pot-type 
combustor  described  earlier.  Additionally,  the  model  will  permit  the 
thermodynamic  characteristics  of  the  present  system  to  be  predicted 
at  operating  conditions  not  examined  experimentally. 

The  specific  objectives  of  the  present  study  can  be  summarized 


as  follows: 


13 


1.  Measure  bath  densities,  bath  compositions,  and  thermal 
performance  for  temperatures  typical  of  metal  combustor 
operation. 

2.  Develop  a  thermodynamic  model  for  the  bath  liquids  and 
compare  the  model  predictions  with  the  measurements. 

3.  Examine  continuous  operation  of  a  complete  steady 
combustor  system  with  varying  thermal  load  requirements 
for  extended  periods  of  time. 

The  study  was  limited  to  the  lithium-sulfur  hexafluoride  reactant 
combination. 
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CHAPTER  II 


THEORETICAL  CONSIDERATIONS 


2.1  Existing  Theory 

Several  techniques  exist  for  the  calculation  of  thermodynamic 
properties  of  multicomponent  solutions.  In  general,  all  techniques 
provide  an  equation  to  represent  the  excess  Gibbs  free  energy  of 

g 

mixing,  f  .  An  excess  thermodynamic  quantity  of  mixing  is  defined  as 
the  deviation  of  that  quantity  from  the  value  given  by  Raoult's  law 
which  governs  an  ideal  mixing  process.  The  excess  Gibbs  free  energy 
of  mixing  can  therefore  be  defined  as: 
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The  two  important  thermodynamic  properties  of  the  mixture, 
enthalpy  and  component  activities,  can  be  obtained  from  the  excess 
Gibbs  free  energy  of  mixing  expression.  The  activity  coefficients  of 
each  component  are  given  by  the  following  exact  thermodynamic  relation: 
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Similarily,  an  expression  can  be  derived  for  the  mixture  enthalpy  as 
a  function  of  the  excess  Gibbs  free  energy  of  mixing. 

The  excess  Gibbs  free  energy  of  mixing  can  be  related  to  the 
excess  entropy  and  enthalpy  of  mixing  as  follows: 
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The  simplest  model  Is  the  ideal  solution  treatment  which  involves 

g 

setting  f  equal  to  zero.  It  is  noted  that  for  an  ideal  solution 
Equation  (2.2)  gives  an  activity  coefficient  of  unity  for  all 
components.  This  approach  was  not  considered  in  the  present  study  since 
the  immiscibility  of  the  system  under  study  indicates  that  it  is 
strongly  nonideal.  Methods  appropriate  to  nonideal  solutions  include 

g 

regular  solution  models  (where  s  is  equated  to  zero) ,  and  athermal 

g 

models  (where  h  is  equated  to  zero) .  The  equations  which  form  the 
thermodynamic  models  contain  empirical  coefficients,  the  number  of 
which  is  dependent  on  the  particular  model  under  consideration  and  the 
number  of  components  in  the  multicomponent  system.  The  advantage  of  the 
models  is  that  the  empirical  coefficients  are  determined  from 
solubility  data  on  the  constituent  binary  systems,  and  in  this  manner 
permit  the  binary  system  data  to  be  extended  to  predict  the  properties 
of  more  complex  systems.  The  thermodynamic  models  available  in  the 
literature  for  this  purpose  are  reviewed  in  the  following  paragraphs, 
Wohl  [27]  provided  a  review  of  three  two-parameter  equations 
referred  to  as  the  -an  Laar,  the  Margules,  and  the  Scatchard-Hamer 
types.  The  van  Laar  equation  has  been  the  most  widely  used,  and  is 
commonly  known  as  the  regular  solution  treatment, 

Wilson  [28]  developed  a  new  type  of  two-parameter  equation 
based  on  the  athermal  model  for  the  excess  free  energy  and 
semitheoretical  considerations.  Applications  of  this  equation  to 
multicomponent  mixtures  were  presented  by  Orye  and  Prausnitz  [29]. 

A  great  disadvantage  of  the  two-parameter  Wilson  equation  is  that  it 
cannot  predict  the  phase  separation  of  immiscible  mixtures,  but  this 
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inadequacy  can  be  removed  by  the  inclusion  of  a  third  empirical 
parameter  forming  a  three-parameter  equation.  The  three-parameter 
Wilson  equation  was  derived  by  Hiranuma  [30]  in  a  manner  which  gave 
some  physical  significance  to  the  third  parameter.  Nagata  and  Yamada 
[31]  modified  the  two-parameter  Wilson  equation  by  expressing  each 
empirical  parameter  as  a  quadratic  function  of  temperature  and  yielding 
a  six-parameter  equation.  Since  the  system  under  study  is  known  to 
be  immiscible,  the  two-parameter  Wilson  equation  was  not  considered 
for  use. 

Renon  and  Prausnitz  [32]  developed  a  three-parameter  equation 
based  on  the  two-liquid  theory  of  binary  mixtures  which  includes 
the  effects  of  the  nonrandoraness  of  mixing.  The  equation  is  referred 
to  as  the  NRTL  (nonrandom,  two-liquid)  equation.  Marina  and  Tassios 
[33,  34]  reduced  the  NRTL  equation  to  a  two-parameter  equation  by 
fixing  the  value  of  negative  unity  for  one  of  the  basic  NRTL  parameters. 
The  resulting  equation  is  referred  to  as  the  LEMF  (local  effective 
mole  fraction)  equation. 

Bruin  [35]  presented  three  additional  two-parameter  equations 
which  were  developed  by  considering  various  approximations  to  the 
excess  Gibbs  free  energy  expression,  Equation  (2.3),  and  employing  basic 
van  Laar  and  Wilson  solubility  parameters.  These  equations  are  known 
as  the  EVL  (extended  van  Laar),  the  enthalpic  Wilson,  and  the  Orye 
equations. 

The  above  review  indicates  that  numerous  thermodynamic  models 
exist  for  the  excess  Gibbs  free  energy  of  a  mixture.  The  models 
were  developed  for  mixtures  composed  of  pure  components  showing 
properties  associated  with  normal  liquids  (see  Appendix  1  of 
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Reference  [38]).  The  question  of  whether  the  models  can  be  applied 
to  the  system  under  study  was  partially  answered  by  Pitzer  [36]  who 
examined  various  alkali  halide-metal  systems,  and  concluded  that  such 
systems  display  thermodynamic  properties  similar  to  those  expected 
for  mixtures  of  normal  non-polar  molecular  liquids.  The  possible 
models  for  the  structure  of  the  alkali  halide-metal  mixtures  were 
examined  by  Lumsden  [37],  who  concluded  that  the  most  plausible  model 
in  the  light  of  Pitzer's  results  is  the  existence  of  metal  atoms  and 
salt  molecules.  It  was  concluded  from  these  findings  that  the  models 
discussed  above  could  be  applied  to  the  present  system  of  alkali  metal 
and  molten  salts. 

The  data  available  on  the  individual  binary  systems  composing 
the  ternary  system  under  study  are  limited  and  do  not  justify  the 
use  of  the  more  complex  models  employing  more  than  two  parameters 
per  binary.  Therefore,  only  the  two-parameter  models  were  considered 
for  use  in  the  present  study.  Most  applications  of  the  models  have 
been  limited  to  systems  composed  of  water  and  organic  compounds 
including  hydrocarbons,  since  the  models  were  developed  primarily  for 
use  in  the  analyses  of  industrial  separation  operations.  The  two 
applications  found  in  the  literature  which  are  related  to  the  present 
study  used  forms  of  the  van  Laar  equation.  Nakanishi  [23]  computed 
the  eutectic  point  of  ternary  molten  salt  mixtures  using  the  Margules 
equation,  which  is  effectively  a  one-parameter  van  Laar  equation 
resulting  from  setting  one  of  the  empirical  parameters  equal  to  unity. 
The  results  indicates  that  the  Margules  equation  computed  the 
eutectic  compoisitions  of  the  mixtures  to  within  1-4  mole  %. 
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Avery  [22]  successfully  employed  the  two-parameter  van  Laar  equation 
in  correlating  the  solubility  characteristics  of  binary  sodium-sodium 
chloride  and  lithium-lithium  fluoride  mixtures .  The  van  Laar  equation 
was  selected  for  use  in  the  present  study  because  it  has  been  used 
to  model  similar  systems,  and  its  simple  mathematical  form 
minimizes  the  difficulty  of  evaluating  the  empirical  parameters.  The 
development  of  the  various  expressions  resulting  from  the  van  Laar 
equation  is  discussed  below. 

2.2  The  Thermodynamic  Model 

2.2.1  Assumptions .  The  major  assumptions  of  the  present 
analysis  are  summarized  as  follows: 

1.  The  van  Laar  expression  for  the  excess  Gibbs  free  energy 
of  mixing  is  used  in  the  liquid  phases. 

2.  The  gas  phase  is  assumed  to  be  an  ideal  mixture. 

3.  The  argon  cover  gas  is  assumed  to  be  insoluble  in  the 
liquids . 

4.  The  change  in  activity  of  the  liquid  phase  components 
with  pressure  is  neglected. 

5.  The  reaction  chamber  pressure  is  taken  to  equal  one 
atmosphere  (1.013  bar). 

6.  The  liquid  phase  densities  are  computed  assuming  no 
volume  change  on  mixing. 

7.  In  the  analysis  of  binary  phase  diagram  data,  it  is 
assumed  that  no  solid  solutions  are  formed. 

2.2.2  Gibbs  Free  Energy  of  Mixing.  The  total  free  energy,  F, 
of  a  homogeneous  solution  is  written  as 
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F  -  E  n  f  °  +  RT  E  n  In  x  +  E  n<  [fE]  (2.4) 

i  i1  1  i  1 

g 

where  f  is  the  excess  free  energy  of  mixing.  If  a  standard  state  is 
chosen  as  the  pure  liquid  component  i,  then  the  activity  is  set  equal 
to  unity  in  this  state,  and  we  have  for  any  other  state 

RT  In  a±  =  f±  -  f±°  (2.5) 


By  differentiating  Equation  (2.4)  with  respect  to  n^  to 
yield  f^,  and  rearranging,  it  can  be  shown  that 

3[n„  fE] 

RT  In  a^  ■  RT  In  x^  +  - ^ -  (2.6) 

1  T,  P,  all  n  (J*i) 

where 


n„ 


£  n 
i  1 


(2.7) 


The  activity  coefficient  is  defined  as  follows: 

a. 


Fixi 


(2.8) 


The  effect  of  pressure  on  the  activity  of  all  components  in  the 
reference  state,  was  neglected  since  the  change  of  activity  with 
pressure  is  small  for  liquids  and  the  system  pressure  is  always 
close  to  atmospheric.  By  substitution  of  Equation  (2.8)  into  Equation 


(2.6),  an  exact  equation  for  the  activity  coefficient  results: 

lT  f 


RT  In  Y  - 


3[nT  fE] 


T,  P,  all  0j  (j/i) 


(2.9) 


2 


As  discussed  in  the  previous  section,  the  van  Laar  expression  for  the 
excess  free  energy  of  mixing  is  employed.  This  expression  is  given 
by  Wohl  [27]  and  Lewis  and  Randall  [34],  and  is  written  for  a  binary 
system  composed  of  components  i  and  j  as  follows: 

„  n.n .b .b .A. . ' 


i  ji  ,1  ij 


(2.10) 


nT(n±bi  +  n^) 

Since  only  two-body  molecular  interactions  were  considered  in  the 

development  of  the  binary  expression,  the  expression  can  be  extended 

to  the  multicomponent  case  by  a  summation  procedure  without  affecting 

the  level  of  approximation  as  follows: 

2  ^  ninjbibjAij ' 

,E  1  j  J  J  J 


"T  £  nibi 


(2.11) 


It  is  noted  that  in  Equation  (2.11)  and  in  all  equations  to  follow, 

V  ■  V  "d  v  ■  V  '  °- 

The  parameters  bi  are  interpreted  as  molecular-volume  or 
quasilattice  parameters,  and  can  be  approximated  by  molar 
volumes  for  some  solutions.  The  parameters  A^  '  are  interpreted 
as  interaction  energies  and  are  a  measure  of  the  deviation  from  ideal- 
solution  behavior.  Typically  both  types  of  parameters  are  positive, 
but  the  interaction  energies  can  be  negative  for  some  solutions. 

Equations  (2.11)  result  in  three  and  six-parameter 
representations  of  binary  and  ternary  systems,  respectively.  The 
binary  expressions  are  reduced  to  two-parameter  equations  and  the 
ternary  expression  to  a  five-parameter  equation  by  the  following 
redefinition  of  parameters: 


2.12 
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bl  b2  bl 

B1  =  b7  B2  "  bj  B3  * 


(2.13) 


The  quantities  defined  by  Equations  (2,12)  and  (2,13)  may  be  functions 
of  temperature.  For  most  systems  the  parameter  is  nearly 
temperature  invariant,  especially  if  it  closely  equals  the  ratio  of 
the  molal  volumes;  and  the  parameter  usually  has  a  small 
temperature  dependence.  It  is  noted  that  the  dissymmetry  of  two  of 
the  three  binary  pairs  composing  the  ternary  system  defines  the 
dissymmetry  of  the  third  pair.  In  equation  form  the  dissymmetry 
relationship  is  stated  as 
B_ 


(2.14) 


Substitution  of  the  parameters  defined  by  Equations  (2.12) 
and  (2.13)  into  Equations  (2.10)  and  (2.11)  results  in  the  following 
expressions  for  the  excess  free  energy  of  the  three  binary  systems 
and  the  ternary  system; 


Binary  1-2: 

fE  _  nln2A12 

Op  (B^n^  +  n2) 

Binary  2-3: 

fE  _  n2n3A23 

ar(B2n2  +  n3) 

Binary  1-3: 

fE  _  nln3A13 

"  nT(B3nl  +  n3} 


nl  +  n2 


n2  + 


ni  + 


(2.15) 


(2.16) 


(2.17) 
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Ternary : 


E  =  nln2B2A12  *  n2n3A23  +  nln3A13 
nT(B3n1  +  %2n2  +  n3^ 


n;L  +  n2  +  n3 


(2.18) 


The  set  of  van  Laar  equations  described  above  are  referred  to 
as  the  unsymmetric  van  Laar  treatment  since  the  equations  include 
effects  of  the  different  molecular  sizes  of  the  pure  components. 

The  symmetric  treatment,  Margules  equation,  results  if  all  the 
parameters  are  set  equal  to  unity. 


2.2.3  Activity  Coefficients.  The  van  Laar  equations  for 
the  activity  coefficients  result  by  substitution  of  Equations  (2.10) 
and  (2.11)  into  Equation  (2.9).  The  general  expression  is  written 
as  follows  prior  to  the  substitution  of  the  parameters  defined  by 


Equations  (2.12)  and  (2.13): 

\  E  E<a  ■  -  i 

In  y  - 


k  ^  j~*  *  ik  2  Aij'^  bibjninj 


E  (nib1)' 


(2.19) 


The  two-parameter  representation  of  a  binary  system  composed  of 
components  1  and  2  is  written  as  follows: 

A12  I  X2  | 

ln  Y1  "  RT~  \B1x1  +  x2| 
r2  Bx  RT  +  x2 


(2.20) 


Similar  expressions  result  for  the  remaining  two  binary  systems  and 
the  ternary  system. 


2.2.4  Mixture  Enthalpy.  The  total  enthalpy  of  a  multicomponent 
mixture  is  defined  as 


H  ■  H.  ,  ,  +  AH  , 

ideal  mix 


(2.21) 


where 


ideal 


=  n  h.C 
.  i  i 


(2.22) 


AH  .  =  AF  .  +  TAS 

mix  mix  mix 


(2.23) 


3AF  . 
mix 


9T  |P,  n, 


(2.24) 


AF  =  F  -  F 
mix  ideal 


(2.25) 


Substitution  of  Equation  (2.4)  into  Equation  (2.25)  results  in 


£  (RTn^  In  x^  +  n^fE) 


(2.26) 


Substitution  of  Equation  (2.26)  into  Equation  (2.24)  results 


-  £  (Rni  In  x±  -  n± 


(2.27) 


Substitution  of  Equations  (2.26)  and  (2.27)  into  Equation 


(2.23)  yields 


AHmix  "  nT  f  "  T  3t~ 


(2.28) 


This  quantity  is  an  excess  quantity  by  definition  since  the  heat 
of  mixing  of  an  ideal  solution  is  zero. 
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The  following  expression  results  for  the  molar  enthalpy  of  the 
solution  by  dividing  Equation  (2.21)  by  n^,  after  substitution  of 
Equation  (2.28): 


h  -  £  x^0  +  fE  -  *  w 


(2.29) 


P,  n. 


The  van  Laar  expression  for  the  excess  free  energy  of  mixing. 

Equation  (2.11),  is  substituted  into  Equation  (2.29)  to  yield 

9b . 


T  EEvjhYij'r 

°  +  — =2 - - - !l  +  T  — -  I 


h  “  E  +  r> 

i  nT  2_  nibi 


P  Vi  J 


E  E* ninj 


-  t  i _ i 


3A  9b  db±l 

bibj  Tt  +  biAij  9T  +  bjAij  9 T 


nT  E  nibi 


(2.30) 


All  partial  derivatives  in  Equation  (2.30)  are  taken  at  constant 
pressure  and  composition. 


2.2.5  Mixture  Density.  The  density  of  the  mixture  was 
calculated  by  neglecting  the  volume  change  of  mixing.  Hildebrand, 
Prausnitz  and  Scott  [39]  presented  an  expression  for  the  change  in 
molal  volume  of  a  component  upon  mixing  that  requires  knowledge  of 

the  internal  pressure,  -gr=-J  ,  of  the  pure  component.  This  correction 

*T 

was  not  used  since  the  internal  pressures  of  the  substances  under 
study  are  not  known;  and  the  calculation  results  in  only  a  small 
correction  to  the  pure  component  value  in  most  systems.  The  density 
of  the  mixture  is  therefore  written  as 


25 


P 


I  MX 

i  1 

I 


(2.31) 


where  and  p.^  are  the  molecular  weight  and  density  of  pure  component 
i,  respectively. 

The  pure  component  properties  used  in  Equations  (2.30)  and 
(2.31)  are  either  taken  from  the  literature  or  estimated  when  the 
necessary  information  was  absent.  Appendix  B  presents  the  sources 
of  the  property  data  and  discusses  the  methods  used  to  estimate 
property  data  not  available  in  the  literature. 


2.2.6  Gas  Phase.  The  gas  phase  is  assumed  to  be  an  ideal 
gas  mixture,  and  therefore  the  fugacity  of  each  species  is  equal 
to  its  partial  pressure,  p^.  The  standard  state  of  a  pure  liquid 
component  is  defined  as  discussed  in  Section  2.2.2.  In  a  state  of 
equilibrium  the  fugacity  of  a  component  is  the  same  in  all  phases. 
Therefore 


where  p^°  is  the  vapor  pressure  of  pure  component  i.  Equation 
(2.19)  is  employed  to  calculate  the  activity  of  component  i  in  the 
solution.  Equations  (2.32)  res-1  K  in  one  equation  for  each  of  the  three 
components,  lithium,  lithium  fluoride  and  lithium  sulfide,  in  the 
ternary  system  under  study. 

The  vapor  pressures  of  the  pure  components  in  the  above  equations 


are  found  from  the  equilibrium  equations  for  vaporization  of  the  pure 
liquids.  Again  a  set  of  three  equations  results: 
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Li (1)  J  Li(g) ;  LiF(l)  t  LiF(g);  Li2S(l)  £  Li2S(g)  (2.33) 

Data  from  the  JANAF  Thermochemical  Tables  [10]  and  the 
experimental  data  of  Berkowitz  and  Chupka  [40]  indicate  the  formation 
of  polymetic  species  in  the  gas  phase.  The  following  equations  are 
written  to  include  the  polymeric  species: 

2Li(g)  £  Li2(g);  2LiF(g)  £  Li2F2(g);  3LiF(g)  £  Li^g)  (2.34) 

The  processes  described  by  Equations  (2.33)  and  (2.34)  can  be 
written  in  general  form  as 


Vi'  Ni  *  Vj"  Nj  (2.35) 

The  general  form  of  the  equilibrium  constant  for  a  process 
described  by  Equation  (2.35)  is  as  follows: 


K. 


(2.36) 


The  equilibrium  constants  of  Equations  (2.36)  can  be  also 
expressed  as 


Af±  -  Af±°  *  RT  In  K±  (2.37) 

At  equilibrium  Af^  ■  0  and  therefore: 

Af^  *  -  RT  In  Kt  (2.38) 

For  the  process  represented  by  Equations  (2.35),  the  following 
expressions  result: 

Kj  -  exp  [(V1'  f±°  -  Vj"  f j°)/RT]  (2.39) 

The  pure  component  properties  used  in  Equations  (2.30)  are  either 
taken  from  the  literature  or  estimated  as  discussed  in  Appendix  B. 


/  -  ar  *.*^  ♦*•*! --.,4  .  ■ 
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The  change  in  the  activities  of  the  pure  liquid  components 
with  pressure  was  neglected,  and  the  activities  were  taken  to  equal 
unity.  The  standard  state  for  the  gaseous  species  was  taken  as 
the  state  where  the  fugacity  equalled  one  atmosphere  (1.013  bar)  to 
permit  the  use  of  tabulated  thermodynamic  data  employing  the  same 
standard  state  [10].  Therefore,  under  the  ideal  gas  assumption  it 
follows  that 

a^g)  “  Pi  (atm)  (2.-40) 

The  partial  pressures  in  Equation  (2.40)  are  denoted  by  p^°  when 
the  equation  is  applied  to  Equations  (2.33)  representing  the 
vaporization  of  the  pure  liquids. 

Substitution  of  Equations  (2.40)  and  (2.32)  into  Equations 
(2.36)  results  in  the  following  expressions  for  the  vapor  pressures 
of  the  various  species  in  units  of  atmospheres: 


'Li  "  ‘S.i  ^i 

Phi  =  K_  ^  a  ^ 

L12  *Li  aLi 

'LiF  “  ^iF  ^iF 

Pl12F2  =  KLi2F2  ^LiF  ^iF 

'Li2S  =  3Li2S 

PLi3F3  “  KLi3F3  ^LiF  ^iF 

(2.41) 

Where  the  activities  are  computed  from  Equations  (2.19)  and  the 
equilibrium  constants  are  computed  from  Equations  (2.39). 

In  addition  to  the  gaseous  species  treated  above,  argon  exists 


in  the  system  as  a  cover  gas.  Thormeier  [41]  investigated  the 
solubility  of  helium  and  argon  in  liquid  sodium,  and  provided  a 
correlation  based  on  a  solubility  parameter  which  permitted 
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prediction  of  the  solubility  of  noble  gases  in  other  liquid  metals. 
After  the  correlation  was  substantiated  by  testing  it  against  the 
data  of  Slotnick,  et  al. ,  [42]  for  the  solubility  of  helium  in 
lithium;  it  was  used  to  predict  that  the  solubility  of  argon  in 
lithium  was  in  the  range  of  10  ^  to  10  ^  on  a  mole  fraction  basis 
at  1200°K.  The  mole  fraction  of  argon  in  a  mixture  of  lithium, 
sodium  and  potassium  fluorides  (46.5-11,5-42.0  mole  percent)  was 
extrapolated  from  data  presented  by  Blander  [43]  to  be  in  the  range 
of  10  ^  to  10'7  at  1200°K.  Based  on  these  results,  the  solubility 
of  the  argon  cover  gas  in  the  liquid  mixture  was  neglected.  For 
the  case  where  the  combustor  is  vented  to  atmospheric  pressure, 
the  total  pressure  of  the  gas  phase  equals  one  atmosphere,  and  the 
partial  pressure  of  the  argon  is  given  by  the  following: 

PA  “  1  *  PLi  "  PLi2  “  PLiF  ‘  PLi2F2  '  PLi3F3  "  PLi2S  (2‘42 

2.2.7  Summary  of  the  Thermodynamic  Model.  The  van  Laar 
thermodynamic  model  developed  in  the  previous  sections  permits  the 
calculation  of  the  excess  Gibbs  free  energy  of  mixing  of  the  ternary 
system  from  data  on  the  solubility  characteristics  of  the 
constitutent  binary  systems.  The  formulation  permits  the  calculation 
of  oth--r  thermodynamic  properties  such  as  liquid  phase  activities, 
densities,  and  enthalpies,  and  gas  phase  partial  pressures.  The 
model  parameters  can  be  calculated  from  binary  solid-liquid  equilibria 
data,  vapor-liquid  equilibria  data,  or  data  on  two  equilibrated 
liquid  phases.  The  procedures  utilized  to  determine  the  parameters 
in  the  present  investigation  are  described  in  Appendix  C,  and  the 
binary  results  are  presented  in  Chapter  III.  The  activity  expressions 


i 


are  utilized  in  the  following  section  to  predict  the  liquid  phase 
solubilities  of  a  steady  metal  combustor  system. 


2.3  Prediction  of  Liquid  Phase  Solubilities  of  a  Steady  Metal 
Combustor 

The  design  of  a  steady  metal  combustor  system  is  based  on  the 
continual  removal  of  the  product-rich  phase  or  the  phase  rich  in 
lithium  fluoride  and  lithium  sulfide,  and  the  continual  addition  of 
lithium  to  replenish  the  fuel-rich  phase.  Separation  of  the  combustor 
bath  liquid  into  the  two  immiscible  phases  occurs  if  the  fugacity 
of  each  component  is  identical  in  each  phase.  Since  the  standard 
state  for  each  component  is  the  same  in  each  phase,  the  condition  for 
phase  separation  implies  that  the  activity  of  each  component  is 
Identical  in  each  of  the  two  phases.  The  two  phases  are  indicated  by 
mole  fractions  and  y^,  and  the  requirement  for  phase  separation 
is  stated  in  equation  form  as 

ai(x±)  ■  a^y^  (2.43) 

Since  the  mole  fractions  in  each  phase  must  sum  to  unity,  two 
additional  equations  result: 

£  x  «  1  (2.44) 

i  1 

£  y.  -  1  (2.45) 

i  1 

Under  steady  state  conditions,  the  fuel  and  oxidizer  flow 
rates  are  in  stoichiometric  proportion,  and  the  products  of  the 
reaction  are  removed  in  a  mole  fraction  ratio  controlled  by  the 


stoichiometry  of  the  reaction.  The  stoichiometry  of  a  chemical 
reaction  can  be  expressed  as 


&  V  *  £  V 


(2.46) 


where  represents  the  chemical  symbol  of  species  i,  and 
represents  the  number  of  moles.  For  such  a  reaction  the  mole  fraction 
ratio  of  any  two  product  species  can  be  expressed  as  follows: 


x  v  " 
i  _  j 

x7  -  v." 
J  3 


(2.47) 


For  the  reaction  under  study  the  stoichiometric  equation  is 


8Li  +  SF6  -*•  6LiF  +  Li2S 


(2.48) 


Therefore,  the  mole  fractions  of  lithium  fluoride  and  lithium 
sulfide  in  the  product-rich  phase  are  related  as  follows: 


*U2S 


(2.49) 


Equations  (2.43),  (2.44),  (2.45),  and  (2.49)  provide  a  set  of 
six  equations  to  solve  simultaneously  for  the  six  unknown  mole 
fractions.  The  procedure  used  to  solve  the  set  of  nonlinear 
equations  is  described  in  Appendix  A.  Following  the  solution  for 
the  unknown  mole  fractions,  the  equations  previously  developed  were 
employed  to  compute  other  mixture  properties  such  as  density  and 
enthalpy,  and  the  gas  phase  composition. 


CHAPTER  III 


DETERMINATION  OF  THE  THERMODYNAMIC  MODEL  PARAMETERS 

3.1  Introduction 

The  thermodynamic  model  for  the  ternary  system  presented  in 
Chapter  II  contains  five  empirical  parameters  which  were  determined 
from  data  on  the  constituent  binary  systems.  The  three  systems  are 
lithium-lithium  fluoride,  lithium  fluoride-lithium  sulfide,  and 
lithium-lithium  sulfide.  The  methods  used  to  obtain  and  analyze  the 
data  on  each  binary  pair  are  described  in  the  following  sections. 

3.2  Lithium-Lithium  Fluoride  Binary 

The  solubility  data  for  the  immiscible  liquid  region  of  the 
lithium-lithium  fluoride  system  presented  by  Dworkin,  et  al. ,  [20] 
were  utilized  to  determine  the  two  empirical  parameters  associated 
with  this  binary  pair.  The  components  of  the  ternary  system  were 
labeled  such  that  subscripts  1,  2,  and  3  represented  components 
lithium,  lithium  fluoride  and  lithium  sulfide,  respectively. 
Employing  this  convention,  the  empirical  parameters  representing  the 
lithium-lithium  fluoride  binary  are  indicated  by  A^  and  B^. 

The  procedure  used  to  determine  the  parameters  is  outlined 
in  Section  C. 2  of  Appendix  C.  Specifically,  the  data  taken  from 
Reference  [20]  was  substituted  into  Equations  (C.17)  and  (C.18)  to 
yield  values  of  the  parameters  over  a  range  of  temperatures.  The 
results  are  shown  in  Figure  3,  along  with  a  linear  least  squares  fit 
for  each  parameter.  The  equations  for  the  parameters  are 
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A12  -  121.53  -  0.05730T  (kJ/mole)  (3.1) 

B1  -  2.5846  -  8.132  x  10-4T  (3.2) 

The  estimated  standard  deviations  of  the  A^  and  parameters  about 
the  regression  lines  are  1.36  kJ  and  0.08,  respectively.  It  is  noted 
that  the  interaction  energy  parameter,  A^,  has  a  stronger  temperature 
dependence  than  the  quasilattice  parameter,  B^.  Since  the  B^ 
parameter  is  associated  with  the  sizes  of  the  molecules,  it  is  not 
expected  to  vary  greatly  with  temperature. 

The  parameters  given  by  Equations  (3.1)  and  (3.2)  were 
employed  with  the  binary  van  Laar  model  to  compute  the  solubility 
characteristics  of  the  immiscible  liquid  region  utilizing  a  procedure 
similar  to  that  described  in  Section  2.3.  In  Figure  4  the  results 
are  compared  to  the  data  taken  from  Reference  [20];  the  agreement 
is  good  over  the  complete  temperature  range  including  the  consolute 
point. 


3.3  Lithium-Lithium  Sulfide  Binary 

Mo  data  were  available  in  the  literature  concerning  the 
lithium-lithium  sulfide  binary  system.  Therefore  tests  were  conducted 
in  this  study  to  determine  the  solubility  characteristics  of  the  region 
of  the  phase  diagram  consisting  of  solid  lithium  sulfide  and  a  lithium- 
rich  liquid. 

The  tests  consisted  of  stabilizing  a  bath  containing  the  two 
components  at  the  desired  temperature,  sampling  the  bath,  and 
subjecting  the  samples  to  the  wet-test  chemical  analysis  procedures 
described  in  Appendix  D.  The  bath  was  contained  in  a  type  316 
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stainless  steel  vessel,  6.27  cm  in  diameter  and  26.0  cm  in  length. 

The  vessel  was  installed  in  a  Lindberg  Hevi-Duty  crucible  furnace, 
type  56622,  controlled  by  a  type  59344  temperature  controller. 
Temperatures  were  measured  with  chromel-alumel  thermocouples 
inserted  into  the  bath,  and  recorded  with  either  a  Leeds  and  Northrup 
model  8686  millivolt  potentiometer  or  a  Leeds  and  Northrup  Speedomax-H 
null-balance  recorder. 

The  lithium  and  lithium  sulfide  reagents  used  in  the  tests  were 
produced  by  the  Lithium  Corporation  of  America  and  Research  Organic/ 
Inorganic  Chemical  Corporation,  respectively.  The  purities  were 
listed  as  99.9%  for  the  lithium  with  the  major  contaminants  being 
sodium  (0.007%)  and  silicon  (0.006%),  and  98%  for  the  lithium  sulfide 
with  the  contaminants  not  specified.  Approximately  30.5  gm  of 
lithium  sulfide  and  166  gm  of  lithium  were  added  to  the  vessel  for 
each  of  the  three  test  sequences. 

Bath  saturation  was  insured  by  maintaining  the  components  at  a 
given  temperature  condition  for  approximately  23  hours.  The  bath 
was  stirred  with  an  argon  bubbler  for  the  first  19  hours  of  the 
period  and  maintained  stagnant  under  an  argon  cover  for  the  remaining 
time . 

The  solubility  test  results  are  plotted  as  mole  fraction  lithium 
sulfide  versus  temperature  in  Figure  5.  The  data  were  employed  to 
compute  the  two  empirical  parameters,  A^  and  ,  associated  with 
this  binary  pair. 

The  solid-liquid  equilibria  expression.  Equation  (C.10), 
developed  in  Appendix  C  was  used  to  evaluate  the  empirical  parameters. 
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A13  -  A  +  B  T 


The  parameter  was  taken  to  be  temperature  invariant,  and  A^  was 
assumed  to  be  a  linear  function  of  temperature  as  follows: 

(3.3) 

In  the  range  where  the  lithium  sulfide  (component  3)  mole  fraction  is 
negligible  compared  to  unity,  Equation  (C.10)  yields  a  one-to-one 
relationship  between  the  parameters  A/B^  and  B/B^  at  each  temperature- 
mole  fraction  condition.  The  equation  was  rearranged  into  the  form 

(3.4) 


Y  »  —  4.  B  y 

i  B3  B3  xi 


where  and  Y^  are  constants  determined  from  the  temperature-mole 
fraction  data.  Least  squares  regression  techniques  were  applied  to  the 
set  of  X^  and  Y^  determined  from  the  data  to  obtain  values  for  A/B3 
and  B/B3. 

Data  on  the  melting  temperature,  heat  of  fusion,  and  heat 
capacity  of  lithium  sulfide  are  required  to  evaluate  the  constants 
X^  and  Y^.  The  melting  temperature  is  known  [15],  but  the  heat  of 
fusion  and  heat  capacity  information  was  not  available  in  the 
literature.  Estimates  of  both  these  quantities  were  made  using  the 
procedures  described  in  Appendix  B.  The  difference  between  the 
partial  molal  heat  capacity  of  lithium  sulfide  in  solution  and  the  molal 
heat  capacity  of  solid  lithium  sulfide  was  taken  to  be  4.18  J/mole 
in  the  temperature  range  of  interest.  The  value  was  estimated  from 
the  estimated  heat  capacities  of  the  pure  liquid  and  solid.  The 
heat-capacity  term  in  Equation  (C.10)  is  dominated  by  the  heat-of- 
fusion  term,  and  therefore  the  heat  capacity  estimates  are  not 


critical  to  the  calculation.  The  heat  of  fusion  of  lithium  sulfide 
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was  estimated  to  be  7.24  KJ/mole.  During  the  analysis  of  the 
lithium  fluoride-lithium  sulfide  binary,  evidence  indicated  that  a 
higher  heat  of  fusion  could  be  associated  with  the  compound.  The 
higher  value  of  49.4  KJ/mole  was  also  carried  through  the  analysis. 

The  data  obtained  on  the  lithium-lithium  sulfide  binary  could 
not  be  used  to  firmly  fix  the  value  of  the  quasilattice 
parameter  B^.  Additional  data  in  the  lithium  sulfide-rich  portion 
of  the  phase  diagram  is  required  to  determine  B^.  Information  could 
be  obtained  on  the  parameter  by  making  a  comparison  to  the  lithium- 
lithium  fluoride  binary.  It  was  expected  that  the  parameter  would 
be  somewhat  smaller  than  the  corresponding  value  of  1.61  (at  1200°K) 
determined  for  the  lithium-lithium  fluoride  binary  because  the  molar 
volume  of  lithium  sulfide  is  larger  than  that  of  lithium  fluoride. 

The  effects  of  the  parameter  B3  on  the  ternary  solubility  predictions 
were  examined.  It  was  found  that  the  experimental  data  was  predicted 
most  accurately  with  a  B^  value  which  yielded  an  interaction  energy 
approximately  equal  to  that  of  the  lithium-lithium  fluoride  binary. 

It  is  not  unexpected  that  the  two  systems  have  similar  interaction 
energies  since  it  is  strongly  suspected,  although  not  verified, 
that  the  lithium-lithium  sulfide  system  has  immiscibility 
characteristics  similar  to  the  lithium-lithium  fluoride  system.  The 
parameter  B^  was  determined  by  equating  the  interaction  energy 
parameter  A^3  to  the  corresponding  parameter  in  the  lithium-lithium 
fluoride  binary,  A^»  at  a  temperature  of  1200°K.  The  values  determined 
in  this  manner  were  1.24  for  the  low  heat  of  fusion  case,  and  1.36 
for  the  high  heat  of  fusion  case.  Both  values  are  lower  than  the 
corresponding  value  for  the  lithium-lithium  fluoride  binary. 
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The  parameters  for  the  low  heat  of  fusion  case  are 


A13  =  88.291  -  0.02975T  (kJ/mole)  (3.5) 

B3  -  1.24  (3.6) 

and  for  the  high  heat  of  fusion  case  are 

A13  -  53.287  -  4.00  x  10-4T  (kJ/mole)  (3.7) 

B3  *  1.36  (3.8) 


The  standard  deviation  of  the  variable  Y^  in  Equation  (3.4)  gives  an 
indication  of  the  precision  of  the  parameter  A  3«  For  the  low  heat 
of  fusion  case,  the  estimated  standard  deviation  of  Y^  is  1.39  kJ/mole 
for  values  ranging  from  47.64  to  59.58  kJ/mole.  For  the  high  heat 
of  fusion  case,  the  estimated  standard  deviation  of  Y^  is  1.44  kJ/mole 
for  values  ranging  from  37.14  to  41.62  kJ/mole.  Changes  in  the 
parameter  B3  affect  the  ternary  predictions  of  lithium  solubility 
in  the  products.  A  B3  variation  of  +0.05  results  in  approximately  a 
+5%  variation  in  the  lithium  solubility  prediction.  The  solid-liquid 
equilibria  results  computed  from  Equation  (C.10)  are  shown  in  Figure 
5  for  the  two  cases  examined. 

3.4  Lithium  Fluoride-Lithium  Sulfide  Binary 

The  solubility  characteristics  of  the  lithium  fluoride-lithium 
sulfide  binary  were  experimentally  examined  by  Faeth  [21]  using  a 
thermal  analysis  procedure.  The  procedure  involved  measuring  the 
thermal  effects  in  the  time-temperature  curves  of  prepared  samples  of 
different  compositions  during  cooling  cycles.  The  thermal  effects 
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consisted  of  constant  temperature  periods  or  hold-ups  at  sample 
temperatures  where  phase  transitions  occurred.  During  the  phase 
transition,  the  temperature  of  the  sample  remained  constant  as  the 
energy  loss  from  the  sample  was  balanced  by  the  enthalpy  change 
associated  with  the  phase  transition. 

The  samples  were  contained  in  a  type  316  stainless  steel 
capsule,  2.22  cm  in  diameter  and  15.2  cm  in  length.  The  capsules 
were  heated  in  a  furnace  to  a  temperature  of  1303 °K  and  allowed  to 
stand  for  a  period  of  five  hours.  The  furnace  was  then  turned  off, 
and  the  time-temperature  curve  of  the  sample  was  recorded.  Two 
thermal  effects  were  observed;  the  effects  represented  the  initial 
solidification  of  lithium  fluoride  from  the  mixture  and  the 
solidification  of  lithium  fluoride  and  lithium  sulfide  at  the 
eutectic  temperature.  Two  tests  were  run  for  each  condition,  and  the 
reproducibility  was  reported  to  be  within  1°K. 

The  data  obtained  from  the  tests  are  plotted  in  Figure  6  as 
mole  fraction  lithium  sulfide  versus  temperature.  For  lithium 
sulfide  concentrations  greater  than  the  eutectic  composition,  only 
one  of  the  two  expected  thermal  effects  was  observed.  The  absence 
of  the  higher  temperature  thermal  effect  associated  with  the 
solidification  of  lithium  sulfide  was  believed  to  be  caused  by 
experimental  difficulties  due  to  supercooling  of  lithium  sulfide  [21]. 

Equation  (C.10)  was  used  in  conjunction  with  the  data  to 
determine  the  empirical  constants  associated  with  this  binary  system, 
and  1*2.  The  required  thermodynamic  data  needed  in  the  analysis 
are  available  in  the  literature  for  lithium  fluoride  and  the  sources 
are  presented  in  Appendix  B.  The  results  yielded  a  value  of  A22 
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equal  to  zero,  indicating  that  the  solution  exhibits  ideal 
solution  behavior  in  the  temperature  and  composition  range  of  the 
data.  With  and  known,  the  value  of  B2  was  determined  by  the 
dissymmetry  relationship.  Equation  (2.14). 

The  solubility  data  also  yielded  information  concerning  the 
heat  of  fusion  of  lithium  sulfide.  Equation  (C.10)  was  used  to  compute 
the  freezing  point  depression  line  extending  from  the  melting  point 
of  pure  lithium  sulfide  to  the  eutectic  point.  The  results  are  shown 
in  Figure  6.  The  difference  between  the  partial  molal  heat  capacity 
of  lithium  sulfide  in  solution  and  solid  lithium  sulfide  was  estimated 
to  equal  an  average  value  of  1.09  J/mole  in  the  temperature  range  of 
interest.  The  freezing  point  depression  line  computed  using  the 
heat  of  fusion  estimated  in  Appendix  B,  7.24  kJ/mole,  does  not  yield 
the  correct  eutectic  point;  the  value  of  the  heat  of  fusion  needed 
to  correctly  compute  the  eutectic  point  was  found  to  equal  49.4  kJ/mole. 
Both  heat  of  fusion  values  were  considered  in  determining  the 
thermodynamic  model  empirical  parameters  as  described  in  Section  3.3 
and  below. 

The  value  of  the  interaction  energy  parameter  is  given  by 

■  0  (kJ/mole)  (3.9) 

It  was  estimated  that  the  value  of  A^  is  accurate  to  +1  kJ/mole. 

The  quaaillatice  parameter  is  given  for  the  low  heat  of  fusion  case 
by 

B2  -  1. 0/(1. 90  -  5.98  x  10_4T) 


(3.10) 
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and  for  the  high  heat  of  fusion  case  by 

B2  -  1. 0/(2. 08  -  6.56  x  10-4T)  (3.11) 

The  accuracy  of  B2  is  fixed  by  the  accuracy  with  which  B^  and  B^  were 
determined  since  it  is  fixed  by  the  dissymmetry  relationship. 

3.5  Summary 

The  five  empirical  parameters  required  in  the  thermodynamic 
model  developed  in  Chapter  II  were  determined.  The  two  parameters 
representing  the  lithium-lithium  fluoride  binary  pair  were  determined 
most  accurately  since  a  greater  amount  of  data  is  available  on  this 
system.  The  two  parameters  associated  with  the  lithium-lithium 
fluoride  binary  are  particularly  important  since  the  ternary  system 
being  considered  in  the  present  study  consists  largely  of  lithium 
and  lithium  fluoride  in  the  composition  and  temperature  range  of 
interest.  Solubility  data  for  the  lithium-lithium  sulfide  and  lithium 
fluoride-lithium  sulfide  binaries  were  less  abundant  because  of  the 
experimental  difficulties  associated  with  the  handling  of  lithium 
sulfide  at  high  temperatures.  Data  generated  for  these  systems  in 
the  present  study  and  by  Faeth  [21]  were  used  to  evaluate  the 
remaining  three  empirical  parameters.  Two  sets  of  parameters  were 
formed,  one  for  each  of  the  two  values  for  the  heat  of  fusion  of  lithium 
sulfide.  Data  on  the  lithium  fluoride-lithium  sulfide  binary  system 
indicated  that  the  heat  of  fusion  of  lithium  sulfide  is  higher  than 
the  estimated  value.  The  predictions  of  the  thermodyanmic  model  as 
applied  to  a  steadily  operating  liquid  metal  combustor  will  be 
compared  to  experimental  data  in  Chapter  V. 


CHAPTER  IV 


EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

4.1  Introduction 

The  experimental  apparatus  used  in  the  study  were  designed  to 
facilitate  the  extraction  of  bath  liquid  samples;  the  measurement 
of  bath  liquid  densities;  and  the  measurement  of  thermal  performance 
under  conditions  where  the  combustion  process  is  conducted  in  a  steady 
manner.  Two  different  combustors  were  employed  in  the  course  of  the 
investigation.  A  radiative  combustor  was  employed  to  generate  the 
bath  liquid  density  data  and  to  obtain  information  required  for  the 
design  of  the  larger  air-cooled  combustor.  The  air-cooled  combustor 
was  used  to  generate  bath  liquid  solubility  and  thermal  performance 
data. 

Each  combustor  was  designed  in  a  breadboard  style  for 
experimental  convenience  in  that  the  reaction  chamber/heat  exchanger, 
product  collection  tank,  exhaust  filter  and  fuel  tank  were  separate 
components,  and  was  not  configured  for  direct  application  in  a 
submersible.  With  the  exception  of  the  exhaust  gas  filter  which  was 
not  used  with  the  radiative  combustor,  both  combustors  consisted  of 
the  following  major  components:  a)  a  reaction  chamber;  b)  a  prod<- 
collection  tank;  c)  a  fuel  tank;  d)  a  fuel  pump;  e)  an  exhaust  gas 
filter;  and  f)  a  control  system.  General  descriptions  of  the 
components  are  given  in  the  following  sections. 
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4.2  Air-Cooled  Combustor  Apparatus 

The  air-cooled  combustor  apparatus  is  shown  schematically 
in  Figure  7,  and  as  a  photograph  in  Figure  8. 

4.2.1  Reaction  Chamber  and  Heat  Exchanger.  The  reaction 
chamber  and  heat  exchanger  were  fabricated  from  type  316  stainless 
steel.  The  reaction  chamber  had  nominal  dimensions  of  0.28  cm,  48.3 
cm  and  26.0  cm  for  the  wall  thickness,  length  and  diameter, 
respectively.  The  annular  heat  exchanger  was  positioned  around  the 
reaction  chamber  by  six  longitudinal  spacers  to  provide  an  annular 
gap  width  of  approximately  0.28  cm. 

A  sump  was  built  into  the  bottom  of  the  reaction  chamber 
to  contain  the  product-rich  liquid.  The  system  was  designed  to 
contain  the  product-rich  liquid  in  a  sump  to  minimize  thermal 
management  problems  associated  with  maintaining  the  products  in  a 
molten  state,  and  to  take  advantage  of  the  superior  heat  transfer 
qualities  of  the  liquid  metal  by  having  molten  metal  in  contact  with 
the  entire  heat-transfer  surface. 

The  product-rich  liquid  was  removed  from  the  reaction  chamber 
through  a  trap  system  which  consisted  of  two  coaxial  tubes  extending 
from  the  top  and  bottom  of  the  chamber  as  shown  in  Figure  7.  The 
shorter  inner  tube  provided  the  path  for  liquid  outflow  to  the  product 
collection  tank.  The  outer  tube  equalized  the  gas  pressures  in  the 
trap  and  the  chamber  ullage  space  by  allowing  argon  cover  gas  to 
flow  from  the  ullage,  through  the  cross-over  port,  and  out  the  trap. 
The  trap  maintained  product  and  fuel  levels  in  the  reaction 
chamber  by  a  simple  hydrostatic  balance.  If  the  top  of  the  product 
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Figure  7  Schematic  of  Air-Cooled  Steady  Combustor  Apparatus 
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layer  is  chosen  as  a  datum  point,  then  the  hydrostatic  pressures  due 
to  the  product  height  in  the  trap  above  this  point  and  the  fuel  height 
in  the  reaction  chamber  above  this  point  are  equal.  The  large 
density  difference  between  the  two  liquids  allows  a  small  product 
height  to  balance  a  much  larger  fuel  height.  The  outer  tube  had  an 
inside  diameter  of  4.75  cm  and  a  wall  thickness  of  0.165  cm,  and  the 
inner  tube  had  an  inside  diameter  of  1.57  cm  and  a  wall  thickness 
of  0.165  cm. 

The  crossover  port  was  positioned  in  a  sump  built  into  the  top 
of  the  reaction  chamber.  The  purpose  of  the  sump  was  to  reduce  the 
heat  losses  from  the  crossover  port  region,  and  to  provide  sufficient 
surface  area  to  mount  an  electrical  guard  heater  used  to  balance 
the  parasitic  heat  loss  from  the  region.  If  the  temperature  of  the 
crossover-port  region  was  permitted  to  drop  below  that  of  the  bath, 
lithium  and  product  vapor  condensation  would  occur  which  could 
over  a  period  of  time  result  in  the  buildup  of  sufficient  liquid  to 
interrupt  the  argon  flow  through  the  trap.  The  sump  temperature  was 
controlled  by  a  1.5  kW  bead  heater  mounted  on  the  walls.  The  heater 
was  manufactured  by  the  Cole-Parmer  Instrument  Company  (model  3116-40) 
and  was  controlled  by  a  Honeywell  temperature  controller  (model 
R7350A) .  The  control  signal  to  the  controller  was  the  output  of  a 
chromel-alumel  thermocouple  which  was  spot-welded  to  the  wall  of 
the  sump. 

A  sampling  port  which  provided  access  to  the  reaction  chamber 
was  mounted  on  the  sump.  The  port  was  primarily  used  to  take  fuel-rich 
liquid  samples  by  lowering  a  sample  capsule  on  a  specially  designed 
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probe  through  the  port,  but  it  also  functioned  as  a  general  inspection 
port  and  a  manual  fuel  addition  port  (for  solid  fuel  sticks). 

Fuel  entered  the  reaction  chamber  through  a  fitting  welded 
into  the  top  of  the  chamber.  A.  low  purge  flow  of  argon  was  provided 
to  the  fitting  to  prevent  diffusion  of  product  vapors  into  the 
fuel  line  during  periods  when  the  fuel  pump  was  inoperative. 

Oxidizer  was  injected  into  the  molten  fuel  through  a  coaxial 
injector  shown  schematically  in  Figure  9.  In  this  type  of  injector 
the  oxidizer  flows  in  an  inner  core  tube  fabricated  from  nickel 
200  alloy  and  is  prevented  from  contacting  the  outer  injector  parts 
fabricated  from  type  316  stainless  steel  by  a  flow  of  inert  argon  gas 
in  the  outer  shroud  tube.  The  injector  was  installed  in  the  reaction 
chamber  by  a  threaded  connection  in  a  boss  mounted  in  the  wall  a 
distance  of  20.3  cm  from  the  top. 

The  cooling  air  entered  the  heat  exchanger  inlet  header  through 
two  inlet  lines  and  was  exhausted  from  the  exit  header  through  two 
exhaust  lines,  with  both  pairs  of  lines  positioned  180c  apart.  The 
double  inlet  and  exit  configuration  was  used  to  ensure  uniform  air 
flow  in  the  heat  exchanger.  The  heat  exchanger  effectiveness  was 
approximately  88%  at  a  bath  temperature  of  1200°K. 

The  exterior  surface  of  the  heat  exchanger  was  equipped  with 
two  electrical  bead  heaters  used  to  melt  the  fuel  prior  to  the  start 
of  the  combustion  process.  The  heaters  provided  an  output  of  1.5  kW 
each,  and  were  manufactured  by  the  Cole-Parmer  Instrument  Company 
(model  3116-40). 

The  reaction  chamber/heat  exchanger  was  connected  to  the 
collection  tank  by  a  product  exhaust  duct.  The  necked-down  section 


Figure  9  Sketch  of  the  Coaxial  Injector 
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of  the  duct  was  designed  to  reduce  the  view  factor  for  thermal 
radiation  from  the  tip  of  the  trap  tube  to  the  cold  collection  tank. 

The  larger  diameter  section  was  designed  to  reduce  the  possibility 
of  product  liquid  contacting  and  freezing  to  the  duct  walls  at  the 
entrance  to  the  product  tank  where  the  walls  are  unavoidably  cooler. 

The  product  sump  and  product  exhaust  duct  wall  temperatures  were 
controlled  by  electrical  bead  heaters  operated  through  a  temperature 
controller.  The  heaters  were  required  tc  balance  the  significant 
heat  loss  through  the  sump/exhaust  region  both  by  conduction  and 
radiation  to  the  collection  tank.  The  heaters  were  manufactured  by 
the  Cole-Parmer  Instrument  Company,  with  a  1.5  kW  heater  (model  3116-40) 
used  on  the  sump  and  necked-down  portion  of  the  exhaust,  and  a  0.4 
kW  heater  (model  3116-10)  used  on  the  larger  diameter  section  of  the 
exhaust.  The  sump  heater  was  operated  by  a  Honeywell  model  R7350A 
temperature  controller.  The  control  signal  to  the  controller  was  the 
output  of  a  chrome 1-alumel  thermocouple  which  was  spot-welded  to  the 
wall  of  the  sump . 

The  exterior  surfaces  of  the  heat  exchanger,  sump,  and  product 
exhaust  duct  were  insulated  with  approximately  a  5  cm  layer  of 
Fiberfrax  Lo-Con  ceramic  fiber  blanket  manufactured  by  the 
Carborundum  Company. 

4.2.2  Fuel  Pump  and  Fuel  Tank.  All  parts  of  the  fuel  pump  and 
fuel  tank  were  fabricated  from  type  316  stainless  steel.  The  fuel 
pump  was  fabricated  from  standard  tubing  with  an  outside  diameter  of 
0.635  cm  and  an  inside  diameter  of  0.457  cm.  The  nominal  dimensions 
of  the  fuel  tank  were  61  cm  in  length  and  27  cm  in  diameter. 
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The  fuel  pump  is  shown  schematically  in  Figure  7.  Argon  gas 
flowed  into  the  pump  through  the  vertical  inlet  line  to  the  footpiece. 
The  gas  then  flowed  as  a  two-phase  mixture  with  liquid  lithium 
through  the  heated  fuel  line  to  the  reaction  chamber.  The  fuel  flow 
rate  was  a  function  of  the  argon  flow  rate,  lift  height,  and  the  fuel 
head  on  the  footpiece  inlet.  Experimentally  determined  pump 
performance  curves  are  shown  in  Figure  10.  The  fuel  line  temperature 
was  maintained  at  approximately  540°K  by  a  flexible,  fibrous-glass 
insulated  heater  in  order  to  prevent  fuel  solidification  in  the  line. 

The  heater  was  manufactured  by  the  Cole-Parmer  Instrument  Company, 

(model  3105-80) ,  and  the  power  output  was  controlled  by  a  variac 
up  to  a  maximum  of  0.38  kW. 

A  manual  fuel  addition  port  was  provided  at  the  top  of  the  fuel 
tank  to  permit  the  addition  of  solid  fuel  which  was  melted  in  the 
tank  and  maintained  at  approximately  540°K.  Power  to  the  tank  was 
provided  by  three  band  heaters  which  were  controlled  by  variacs . 

The  heaters  were  manufactured  by  the  Chromalox  Corporation,  (model 
HB-1120) ,  and  provided  a  maximum  output  of  2.02  kW  for  each  heater. 

The  exterior  surfaces  of  the  fuel  tank  were  insulated  with  approximately 
a  5  cm  layer  of  Fiberfrax  Lo-Con  blanket. 

4.2.3  Product  Collection  Tank.  The  product  collection  tank 
was  fabricated  from  type  304  stainless  steel.  The  nominal  dimensions 
of  the  tank  were  91  cm  in  length  and  45  cm  in  diameter. 

The  product  exhaust  duct  was  connected  to  the  collection  tank 
lid  with  a  bolted  flange  sealed  by  a  metal  0-ring  (Advanced  Products 
P/N  4750-5-7-SP) .  The  lid  was  the  only  part  of  the  apparatus  bolted 
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Figure  10  Performance  Curves  of  the  Argon  Lift  Pump 
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to  the  frame,  and  functioned  as  a  single-point  mount  for  the  reaction 
chamber/heat  exchanger  and  the  collection  tank.  The  collection  tank 
was  connected  to  the  lid  with  a  bolted  flange  sealed  with  a  Viton  0-ring. 

The  product  collected  tank  was  cooled  by  coils  soldered  to  the 
tank  wall,  using  water  as  the  coolant. 

Two  ports  were  provided  for  access  to  the  collection  tank  as 
shown  in  Figure  7.  The  top  port  was  used  to  take  product-rich  liquid 
samples.  The  sample  was  retracted  and  allowed  to  cool  in  the  sample 
chamber  prior  to  final  removal.  The  lower  inspection  port  was  used 
to  observe  the  product  level  in  the  tank,  and  to  provide  access  to 
the  tank  for  physical  manipulation  of  the  product  as  the  situation 
required. 

4.2.4  Exhaust  Filter.  All  argon  flows  entered  the  collection 
tank  and  then  passed  through  an  exhaust  filter  prior  to  being  exhausted 
into  a  hood.  The  filter  nominal  dimensions  were  16.5  cm  in  diameter 
and  45.5  cm  in  length.  The  filter  element  consisted  of  two  sizes  of 
type  316  stainless  steel  wire  mesh  cloth.  The  inner  wraps  were 

40  x  40  mesh  cloth  followed  by  outer  wraps  of  200  x  200  mesh  cloth. 

No  attempt  was  made  to  optimize  the  filter  design. 

4.2.5  Control  System.  Figure  11  shows  a  schematic  diagram 

of  the  flow  control  system  for  the  apparatus.  Details  of  some  of  the 
instrumentation  shown  in  the  figure  are  discussed  in  the  next  section. 

The  various  purge  flows  shown  in  Figure  11  were  employed  to  prevent 
clogging  of  the  gauge  lines  in  the  system  by  condensed  lithium  vapor. 

The  purge  flows  also  functioned  to  alleviate  air  contamination  of 
the  system  when  sample-port  caps  were  removed.  All  flows  in  the 
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system  were  metered  with  rotameters  with  the  exception  of  the  oxidizer 
and  coolant  air  flows  which  were  metered  with  critical  flow  orifices. 

4.3  Radiative  Combustor  Apparatus 

The  fuel  tank,  fuel  pump  and  control  system  were  common  to 
both  the  radiative  and  air-cooled  combustors  and  are  described  in 
Section  4.2. 


4.3.1  Reaction  Chamber.  The  reaction  chamber  was  fabricated 
from  six  inch  standard  I.P.S.,  schedule  10,  type  316  stainless  steel 
pipe,  nominally  16.5  cm  in  diameter  and  46  cm  in  length.  The 
energy  released  by  the  reaction  was  removed  by  convection  and 
radiation  from  the  side  walls  of  the  reaction  chamber.  The  coaxial 
injector  used  to  inject  oxidizer  into  the  reaction  chamber  and  the 
trap  system  used  to  remove  reaction  products  are  similar  to  the  air¬ 
cooled  combustor  and  a  discussion  of  these  components  is  found  in 
Section  4.2.1.  The  product-rich  liquid  was  contained  within  the 
main  reaction  chamber  at  a  depth  of  15  cm  in  early  tests  and  7.5  cm 
in  later  tests.  The  large  product-rich  liquid  volume  permitted 
accurate  measurements  of  the  density  of  this  phase  to  be  obtained. 

Three  1.5  kW  electrical  bead  heaters,  manufactured  by  the 
Cole-Parmer  Instrument  Company  (model  3116-40) ,  were  installed  on  the 
reaction  chamber  walls  in  order  to  heat  the  fuel  to  the  molten  state 
prior  to  ignition.  Two  of  the  heaters  were  turned  off  following 
ignition,  but  the  heater  installed  in  the  trap  area  was  controlled 
manually  during  a  test  to  balance  the  parasitic  heat  loss  from  this 


region. 
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The  exterior  surfaces  of  the  reaction  chamber  were  insulated 
with  approximately  a  5  cm  layer  of  Fiberfrax  Lo-Con  ceramic  fiber 
blanket.  The  energy  release  rate  from  the  reaction  chamber  was 
controlled  by  manual  adjustment  of  the  amount  of  insulation  on  the 
reaction  chamber  walls. 

4.3.2  Product  Collection  Tank.  The  product  collection  tank 
was  fabricated  from  10  inch  standard  I.P.S.,  schedule  10,  type  316 
stainless  steel  pipe,  nominally  27.3  cm  in  diameter  and  61  cm  in 
length.  The  tank  was  fitted  with  removable  end  plates  on  both  ends 
to  facilitate  post-run  removal  of  the  reaction  products.  The 
reaction  chamber  was  mounted  to  the  collection  tank  by  a  bolted 
flange  sealed  with  a  ceramic  fiber  gasket. 

4.4  Ins  trumentation 

4.4.1  Flow  Measurement.  All  argon  flows  to  the  system  were 
metered  with  rotameters  manufactured  by  Matheson  Gas  Products. 

Rotameter  tube  sizes  varied  depending  on  the  magnitude  of  the  particular 
flow  being  measured.  The  largest  argon  flow  to  the  system  was  through 
the  coaxial  injector  where  a  No,  604  tube  was  utilized.  All 
rotameters  were  calibrated  with  a  Precision  Scientific  Company  wet-test 
meter  (0.10  cu.  ft. /rev.). 

The  oxidizer  flow  was  metered  and  controlled  with  a  critical- 
flow  orifice-pressure  regulator  combination.  The  pressure 
regulator  was  a  Matheson  model  9  two-stage  regulator  with  a  0-300  psig 
output  capability.  The  pressure  upstream  of  the  critical  flow  orifice 
was  measured  with  a  Heise  absolute  pressure  gauge  with  a  0-300  psia 
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capacity.  The  orifice  was  calibrated  with  a  Precision  Scientific 
Company  wet-test  meter  (0.25  cu.  ft. /rev.) 

The  coolant  air  flow  measurement  system  was  divided  into  low 
and  high  range  sections.  The  low-range  section  consisted  of  a 
calibrated  1.78  cm  (0.70  inch)  diameter  flow  orifice  manufactured 
by  the  American  Meter  Company  (model  22092) .  The  pressure  drop 
across  the  orifice  was  measured  with  a  Heise  absolute  pressure  gauge 
with  a  0-300  psia  capacity.  The  high-range  section  consisted  of 
two  critical  flow  orifices  installed  in  the  heat  exchanger  inlet 
lines.  The  low-range  and  high-range  systems  were  designed  so  that 
choked  flow  occurred  in  the  orifices  when  the  air  flow  exceeded  the 
capacity  of  the  low-range  section.  Flow  through  the  critical  flow 
orifices  was  then  metered  by  measuring  the  upstream  pressure  on  the 
line  pressure  gauge  described  earlier.  Flow  through  the  system  was 
controlled  with  a  pressure  regulator  manufactured  by  the  Airline 
Hydraulics  Corporation  (Componetrol  model  420324-5) .  Both  the  low- 
range  and  high-range  sections  were  calibrated  with  a  positive- 
displacement  rotary  meter  manufactured  by  the  Roots-Connersville 
Blower  Corporation  (0.2633  cu.  ft. /rev.). 

A  4.44  cm  diameter  flow  orifice  was  installed  in  the  exit 
plane  of  each  air  exhaust  line.  The  pressure  drop  across  the 
orifices  in  each  line  was  measured  to  allow  corrections  to  be  made 
should  the  flows  in  each  leg  become  unmatched  during  the  thermal 
performance  tests.  Pressures  were  measured  with  an  Ellison  4-inch 
inclined  manometer.  The  total  pressure  drop  through  the  heat 
exchanger  and  exhaust  lines  were  measured  with  a  Merriam  30-inch 
manometer  (model  RC  1418) . 
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4.4.2  Temperature  Measurement.  Temperatures  were  measured 
with  chrome 1-alumel  thermocouples.  The  reaction  chamber  bath 
temperature  was  measured  at  two  different  depths  with  four  grounded, 
type  316  stainless  steel-sheathed  thermocouples  manufactured  by  the 
Conax  Corporation  (P/N's  K-3165512-G-T1-12  and  K-3165512-G-T1-18) . 

The  exhaust  air  temperatures  were  also  measured  with  Conax  grounded, 
type  316  stainless  steel-sheathed  thermocouples  (P/N  K-316556-G-PJ- 
MX-062-A-3).  Two  thermocouples  were  positioned  in  the  center  of  each 
exhaust  leg.  The  surface  temperature  of  various  components  in  the 
apparatus  were  measured  by  spot-welded  thermocouples  fabricated 

from  wire  manufactured  by  the  Thermo-electric  Company,  type  G/G-26-K. 
All  thermocouples  were  manufactured  to  the  American  National  Standards 
Institute  limits  of  error.  In  the  temperature  range  of  use,  the  error 
limit  is  +0.75  percent. 

Temperatures  were  recorded  every  eighteen  seconds  on  three 
Leeds  and  Northrup  Speedomax-H  twelve-point,  null-balance  recorders. 
The  extension  leads  and  recorders  were  calibrated  with  a  Leeds  and 
Northrup  model  8686  millivolt  potentiometer.  The  temperature- 
millivolt  conversion  values  were  taken  from  National  Bureau  of 
Standards  tables. 

4.4.3  Fuel  Level  Measurement.  The  fuel  levels  in  both  the 
reaction  chamber  and  fuel  tank  were  measured  with  argon  bubblers 
located  at  known  positions  in  the  respective  vessels.  A  low  flow  of 
argon  (approximately  0.03  kg/hr)  was  supplied  to  the  bubbler, 

and  the  inlet  pressure  was  measured  on  a  manometer.  The  frictional 
losses  in  the  bubbler  were  negligible  and  therefore  the  manometer 
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measured  the  fuel  head  at  the  bubbler  exit.  An  Ellison  8-inch 
manometer  and  a  Merriam  35-inch  manometer  (model  RB  1418-1)  were 
used  to  measure  bubbler  inlet  pressures  in  the  reaction  chamber  and 
fuel  tank,  respectively.  Knowledge  of  the  bath  liquid  densities 
allowed  the  pressure  measurements  to  be  converted  to  depth  readings. 

4.4.4  Sampling  System.  Fuel-rich  bath  liquid  samples  were 
obtained  for  analysis  by  lowering  a  sample  cup  held  by  a  long-handled 
sampling  probe  through  the  fuel  sampling  port  at  the  top  of  the 
apparatus.  The  probe  was  fabricated  from  type  316  stainless  steel 
and  consisted  of  a  support  for  the  sample  cup  and  a  movable  ball 
which  acted  as  a  filling  valve.  The  sample  cup  was  fabricated  from 
nickel  200  alloy,  and  had  an  outside  diameter  of  1.3  cm  and  a  length 
of  2.5  cm.  During  the  course  of  the  investigation  it  was  found  that 
sample  cups  fabricated  from  either  stainless  steel  or  low-carbon 
steel  resulted  in  sample  contamination  believed  to  be  caused  by 
surface  corrosion  of  the  material  by  the  lithium  sulfide.  The  probe 
and  sample  cup  apparatus  were  similar  to  that  used  by  Avery  [22]. 

The  sample  cup  was  placed  in  the  sampling  probe  and  sealed 
with  the  ball  to  prevent  contamination  prior  to  being  filled  with  the 
melt.  After  the  sample  was  taken,  the  cup  was  again  sealed  with  the 
ball  to  prevent  sample  contamination  during  retrieval. 

Product-rich  bath  liquid  samples  were  obtained  by  sampling  the 
product  flowing  out  the  trap.  The  sample  was  taken  by  extending  a 
sampling  cup  through  the  sample  chamber  into  the  product  flow  path. 
Figure  3  shows  the  product  sampling  cup  in  sampling  position.  The  cup 
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was  fabricated  from  type  316  stainless  steel  and  had  a  semicircular 
cross  section  with  an  outside  diameter  of  3.18  cm  and  a  length  of 
7.62  cm. 


4.5  Experimental  Procedure 

4.5.1  Reactants ■  The  lithium  fuel  used  in  the  experimental 
program  was  produced  by  the  Lithium  Corporation  of  America.  The  fuel 
was  packed  in  two-pound  cans  and  sealed  under  argon.  The  purity  is 
listed  in  99.9%  with  the  major  impurity  being  sodium. 

The  sulfur  hexafluoride  oxidizer  was  produced  by  Matheson  Gas 
Products.  The  gas  was  contained  in  a  cylinder  as  a  liquid  under  its 
own  vapor  pressure  (2200  kPa  at  294°K).  The  gas  was  certified  purity 
grade  with  a  listed  purity  of  99.8%. 

4.5.2  Pre-Test  Preparation  Procedure.  Prior  to  the  start  of 
the  combustion  process,  the  reaction  chamber  including  the  sumps  and 
outflow  duct  were  heated  electrically  to  approximately  1030°K. 
Ignition  could  be  initiated  at  any  temperature  above  the  melting 
point  of  lithium,  453°K,  but  the  higher  temperature  was  chosen  to 
shorten  the  low-temperature  transient  period  of  the  run.  The  fuel 
tank  and  fuel  line  were  heated  electrically  and  stabilized  at  a 
temperature  of  540°K.  During  the  heat-up  period  all  argon  flows 

to  the  system  were  operational  including  a  purge  flow  of  argon 
through  the  core  tube  of  the  coaxial  injector. 

4.5.3  Normal  Operating  Procedure.  The  combustion  process 
was  started  by  actuation  of  a  solenoid  valve  which  supplied  oxidizer 
to  the  injector.  The  argon  purge  flow  to  the  injector  core  tube  was 
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terminated  immediately  after  start.  The  bath  temperature  was  allowed 
to  increase  to  the  desired  operating  value ,  and  the  oxidizer  and 
air  flow  were  then  matched  to  stabilize  the  temperature. 

The  oxidizer  flow  was  adjusted  to  the  desired  power  level, 
and  the  argon  flow  to  the  lift  pump  was  adjusted  to  yield  the 
stoichiometric  fuel  flow  rate.  The  fuel  flow  rate  was  not  measured 
directly,  but  was  monitored  by  observing  the  fuel  level  in  the 
reaction  chamber  through  the  use  of  the  fuel  level  manometer.  Correct 
matching  of  the  fuel  and  oxidizer  flow  rates  was  indicated  by  a 
constant  fuel  level  in  the  reaction  chamber.  The  fuel  level  in  the 
remote  fuel  tank  was  maintained  by  the  manual  addition  of  solid  fuel 
to  the  tank.  The  coolant  air  flow  to  the  system  was  adjusted  at 
each  power  level  to  maintain  the  bath  temperature  at  the  desired 
value. 

The  combustor  was  operated  at  a  particular  power  level-bath 
temperature  condition  for  sufficient  time  to  ensure  the  establishment 
of  thermal  and  chemical  equilibrium  in  the  reaction  chamber.  The 
design  of  the  combustor  facilitated  the  establishment  of  chemical 
equilibrium  since  the  product-rich  liquid  contained  in  the  combustor 
sump  could  be  exchanged  relatively  rapidly  (approximately  0.54  hours 
at  a  power  level  of  10  kW) .  Additionally,  the  fuel-rich  phase  was 
easily  saturated  with  products  since  the  bath  was  stirred  by  the 
oxidizer  jet,  and  the  combustion  process  was  occurring  within  the 
fuel-rich  phase.  Both  fuel-rich  and  product-rich  liquid  samples 
were  obtained  over  a  range  of  temperatures  by  gradually  increasing 
the  bath  temperature  over  the  course  of  the  run. 


The  oxidizer  flow,  shroud  argon  flow,  fuel  pump  argon  flow, 
coolant  air  flow,  fuel  tank  and  reaction  chamber  fuel  levels, 
combustor  pressures,  and  system  temperatures  were  continually 
monitored  and  recorded  throughout  a  test.  The  data  were  later 
reduced  to  determine  system  performance. 

The  combustion  process  was  terminated  by  actuating  the 
oxidizer  solenoid  valve  to  the  closed  position.  An  argon  purge  flow 
through  the  injector  core  tube  was  initiated  just  prior  to  termination 
of  the  oxidizer  flow. 

4.5.4  Post-Test  Shutdown  Procedure.  After  the  oxidizer  flow 
was  terminated,  the  coolant  air  flow  was  reduced  to  a  low  value  to 
allow  gradual  cool-down  of  the  apparatus.  The  product  collection 
tank  was  removed  from  the  apparatus  after  the  system  had  completely 
cooled,  and  the  reaction  products  were  removed  from  the  collection 
tank.  The  product  exhaust  duct  region  of  the  apparatus  was  inspected, 
and  the  collection  tank  was  reinstalled  on  the  apparatus.  The 
injector  and  fuel  level  bubbler  were  removed  and  inspected;  the 
parts  were  replaced  if  necessary. 

4.5.5  Bath  Sampling  Procedure.  Product-rich  and  fuel-rich 
bath  samples  were  extracted  with  the  apparatus  described  in  Section 
4.4.4. 

A  product-rich  sample  was  allowed  to  cool  in  the  sample  chamber 
under  an  argon  purge  prior  to  final  removal  from  the  apparatus.  The 
sample  was  then  cracked  loose  from  the  sample  cup  into  a  stainless 
steel  tray,  and  transferred  to  a  sealed  glass  container. 
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A  fuel-rich  sample  was  removed  from  the  reaction  chamber 
through  the  sample  port  at  the  top  of  the  combustor.  The  combustion 
between  the  lithium  on  the  outside  of  the  sample  cup  and  the  ambient 
air  was  extinguished  and  the  residue  removed  from  the  outside  of 
the  cup  by  quenching  in  water.  The  sample  cup  was  held  in  the 
sample  probe  during  this  period  and  sealed  by  the  probe  ball  valve 
to  prevent  sample  contamination.  The  sample  cup  was  then  removed 
from  the  probe  and  transferred  to  a  sealed  glass  container. 

The  sealed  glass  containers  containing  the  samples  were 
evacuated  in  the  antechamber  of  a  vacuum/glove  box.  The  containers 
were  then  transferred  to  the  main  chamber  of  the  glove  box  and 
stored  under  a  low  humidity,  argon  atmosphere  for  future  chemical 
analysis. 


4.5.6  Bath  Sample  Analysis  Procedure.  The  fuel-rich  samples 
were  prepared  for  analysis  by  cleaning  the  remaining  scale  from  the 
outside  of  the  cups  with  fine  emergy  cloth.  The  polishea  cups  were 
weighed  with  an  analytical  balance  manufactured  by  Ainsworth, 
Incorporated,  type  L  C  28948.  The  cups  were  reweighed  following  the 
analysis  of  the  samples  to  obtain  the  sample  weight  by  difference. 

The  product-rich  samples  were  prepared  for  analysis  by 
grinding  the  entire  sample  to  an  extremely  fine  powder  with  a  model 
7875  Micro-Mill  manufactured  by  The  Chemical  Rubber  Company.  A  portion 
of  the  powdered  sample  was  weighed  on  the  analytical  balance  described 
earlier  and  then  subjected  to  the  analysis  procedure.  All  the 
cleaning,  grinding  and  weighing  operations  were  conducted  within  the 
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low  humidity,  argon  atmosphere  glove  Sjx  to  prevent  sample 
contamination. 

The  details  of  the  wet-test  chemical  analysis  procedures  used 
to  analyze  the  samples  are  described  in  Appendix  D. 

4.5.7  Bath  Density  Measurement  Procedure.  The  densities  of 
both  the  fuel-rich  and  product-rich  liquids  were  obtained  over  a 
range  of  temperatures.  The  radiative  combustor  described  in  Section 
4.3  was  utilized  to  generate  the  density  data.  The  product-rich 
liquid  was  contained  in  the  bottom  of  the  cylindrical  reaction 
chamber  in  the  radiative  combustor.  The  design  permitted  direct 
access  to  both  the  fuel-rich  and  product-rich  liquids  from  the 
sampling  port  provided  in  the  top  of  the  combustor. 

The  measurements  were  made  by  lowering  a  bubbler  (similar  to 
the  fuel  level  bullers  previously  described  in  Section  4.4.3)  through 
the  sampling  port  into  the  bath  with  the  pressure  being  measured 
every  1.27  cm.  The  pressure  was  measured  with  a  Merriam  manometti 
(model  RC  1418)  filled  with  Merriam  unity  oil  (specific  gravity  =  1.00). 
The  reaction  chamber  pressure  was  referenced  to  the  low  pressure 
side  of  the  manometer;  therefore,  the  manometer  measured  liquid  head 
above  the  bubbler  tip.  A  series  of  readings  were  obtained  as  the 
bubbler  was  lowered,  with  approximately  15-22  readings  taken  in  the 
fuel-rich  liquid  and  5-10  readings  in  the  product-rich  liquid. 

A  linear  least  squares  curve  fit  was  made  to  the  pressure  readings; 
the  slope  of  this  curve  yields  the  density.  A  change  in  slope  was 
noted  when  the  product-rich  liquid  was  entered,  and  the  point  where 
the  slope  changed  provided  a  measure  of  the  product-rich  liquid  level. 
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4.5.8  Thermal  Performance  Measurement  Procedure.  The  thermal 
performance  measurements  were  made  by  measuring  the  enthalpy  rise 
and  flow  rate  of  the  coolant  air  as  a  function  of  oxidizer  flow 
rate  for  a  bath  temperature  of  1197+3°K.  The  temperature  was  chosen 
because  it  is  in  the  middle  of  the  normal  operating  range  of  the 
combustor.  The  single  temperature  was  examined  because  the  thermal 
performance  is  not  significantly  affected  by  bath  temperature 
in  the  temperature  range  of  interest. 

For  a  given  oxidizer  flow  rate,  the  coolant  air  flow  rate  was 
adjusted  to  maintain  the  constant  bath  temperature.  The  air  and 
wall  temperature  of  each  exhaust  leg  were  measured  and  recorded. 

The  exhaust  air  temperatures  were  corrected  for  errors  caused  by 
conduction  and  radiation  to  the  walls  by  employing  the  procedures 
presented  by  Moffat  [44].  The  dominant  effect  was  radiation,  with 
corrections  in  the.  range  of  1-5%  resulting  for  the  different  flow 
conditions.  Minor  flow  mismatches  in  the  exhaust  legs  were  accounted 
for  by  applying  corrections  generated  from  data  on  the  pressure  drops 
across  the  exhaust  orifices.  Flow  mismatches  on  the  order  of  5% 
were  encountered  at  the  highest  flows.  The  combustor  was  operated 
for  a  period  of  one  hour  at  each  power  level.  Minor  variations  in 
the  air  flow  rate  and  air  exit  temperature  over  the  hour  period  were 
accounted  for  by  taking  mean  values.  The  enthalpy  change  of  the  air 
was  obtained  by  integrating  the  cubic  expression  for  the  specific 
heat  at  constant  pressure  of  air  given  in  Reference  [45]  from  the 
inlet  temperature  to  the  exit  temperature.  The  enthalpy  change  of 
each  exhaust  leg  was  multiplied  by  the  flow  rate  of  that  leg  and 
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then  added  to  yield  the  total  power  input  to  the  coolant  air  load. 


The  electrical  power  inputs  to  the  system  were  independent  of 

load  and  were  therefore  included  with  the  parasitic  heat  loss  from  the 

system  which  includes  all  heat  flows  not  transferred  to  the  working 

fluid  of  the  heat  exchanger.  Both  the  energy  release  rate  of  the 
« 

combustor  and  the  net  average  parasitic  heat  loss  from  the  combustor 
were  obtained  from  a  plot  of  heat  exchanger  load  in  kW  versus 
oxidizer  flow  rate  in  kg/hr.  The  slope  of  the  curve  yielded  the 
energy  release  rate  in  kW-hr/kg  of  oxidizer.  The  negative  intercept 
on  the  load  axis  yielded  the  net  average  parasitic  heat  loss  from 
the  system  in  kW.  The  gross  parasitic  heat  loss  equals  the  net 
parasitic  heat  loss  plus  the  electrical  input  to  the  system  of 
approximately  1.9  kW. 


CHAPTER  V 


EXPERIMENTAL  AND  THEORETICAL  RESULTS 

5.1  General  Experimental  Observations 

5.1.1  Radiative  Combustor  Tests.  The  objectives  of  the 
radiative  combustor  test  program  were  to  generate  bath  liquid 
density  data  and  to  examine  several  design  details  which  had  to  be 
resolved  before  the  design  of  the  more  complex  air-cooled  combustor 
could  be  completed.  A  total  of  eight  tests  were  conducted,  as 
summarized  in  Table  1. 

Test  LSR-1  indicated  that  a  guard  heater  was  required  in  the 
trap  area  of  the  combustor  to  balance  the  parasitic  heat  loss  from 
this  zone,  and  larger  crossover  ports  were  needed  in  the  trap  shroud 
tube.  Product  outflow  was  not  achieved  during  the  test  because  the 
trap  area  could  not  be  maintained  at  a  temperature  above  the 
solidification  temperature  of  the  product.  These  design  changes 
proved  to  be  satisfactory  during  test  LSR-2,  but  an  additional 
difficulty  was  encountered  with  the  type  321  stainless  steel 
hydroformed  bellows  used  as  the  product  exhaust  duct.  The  bellows  was 
found  to  be  plugged  with  product  material  and  severely  corroded.  To 
eliminate  the  problem,  the  trap  area  was  redesigned  to  provide  a 
cylindrical  product  exhaust  duct  fabricated  from  type  316  stainless 
steel,  along  with  a  modified  inner  trap  tube  which  included  a  lip 
to  direct  the  product  flow  away  from  the  exhaust  duct  walls.  No 
problems  were  encountered  with  the  modified  combustor  during  test 
LSR-3  which  was  run  for  a  total  of  7.1  hours.  The  final  design 


Table 


LSR-8  4.5  0.0-2,77  1252  Demonstrated  extended  high  power  operation.  Bath  densities 

_  measured  at  1177°K. 


70 


modifications  to  the  combustor  were  made  following  test  LSR-6. 

Minor  difficulties  were  encountered  with  product  liquid  freezing  to 
the  exhaust  duct  walls  in  the  region  where  the  duct  was  connected  to 
the  product  collection  tank.  The  duct  diameter  was  increased  in  this 
region  to  minimize  the  possibility  of  product  flow  contacting  the 
walls.  The  inner  trap  tube  was  also  shortened  to  provide  a  lower 
product  level  in  the  combustor. 

The  first  attempt  to  obtain  bath  density  data  was  made  during 
test  LSR-4,  but  a  malfunction  of  the  density  measurement  apparatus 
caused  the  test  to  be  terminated  early.  The  density  measurement 
apparatus  was  modified  and  functioned  properly  during  test  LSR-5 
which  was  run  for  a  total  of  7.1  hours.  The  long-term  operation 
experienced  during  LSR-5  provided  further  substantiation  of  the  design 
modifications  described  earlier.  Additional  density  data  were 
obtained  during  all  subsequent  radiative  combustor  tests. 

The  operation  of  the  argon  lift  pump  used  to  supply  fuel 
to  the  reaction  chamber  was  evaluated  during  the  final  three 
radiative  combustor  tests,  LSR-6,  LSR-7  and  LSR-8.  Prior  to  test 
LSR-6,  the  fuel  had  been  added  to  the  combustor  manually  as  solid 
sticks.  It  was  found  that  the  performance  of  the  pump  was  adversely 
affected  by  small  positive  pressure  differentials  between  the 
reaction  chamber  and  the  fuel  tank,  and  the  argon  exhaust  line  from 
the  product  collection  tank  was  enlarged  to  decrease  the  combustor 
pressure. 

Test  LSR-8  was  conducted  to  demonstrate  long-term  operation 
at  power  levels  greater  than  11  kW.  Figure  12  illustrates  a  portion 
of  the  data  obtained  from  the  test.  A  power  level  greater  than  11  kW 
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Figure  12  Radiative  Combustor  Run  Data  (Test  LSR-8) 
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was  maintained  for  over  three  hours,  with  the  power  exceeding  13  kW 
for  one  hour.  An  idel  period  was  necessitated  midway  through  the 
high  power  period  while  repairs  were  made  to  a  failed  oxidizer 
control  regulator.  The  rather  large  fluctuations  in  the  combustor 
wall  temperature  and  the  product  exhaust  duct  wall  temperature 
shown  in  Figure  12  were  caused  by  the  fact  that  the  combustor 
was  radiatively  cooled  which  resulted  in  the  power  output  of  the 
combustor  being  strongly  affected  by  small  fuel  level  changes.  Also 
shown  in  the  figure  are  the  oxidizer  flow  rate  as  a  function  of 
time  and  the  corresponding  injector  shroud  argon  flow. 

5.1.2  Air-Cooled  Combustor  Tests.  Table  2  summarizes  the 
ten  air-cooled  combustor  tests  which  were  conducted.  The  main 
objectives  of  the  tests  were  to  generate  bath  liquid  solubility  and 
thermal  performance  data,  and  to  demonstrate  long-term  system 
operation  under  steady-state  conditions. 

Tests  PPC-1  and  PPC-2  indicated  the  need  for  a  change  to  the 
combustor  trap  design.  It  was  found  that  the  top  of  the  combustor, 
including  the  region  of  the  shroud  tube  crossover  port,  could  not 
be  maintained  at  bath  temperature;  which  resulted  in  condensed 
lithium  and  product  vapors  forming  a  liquid  slug  inside  the  trap 
shroud  tube.  The  slug  formation  was  enhanced  by  the  exhaust  argon 
flow  causing  condensed  vapors  to  be  directed  into  the  tube,  the  high 
surface  tension  of  liquid  lithium  permitting  the  size  of  the  slug 
to  be  large,  and  the  fact  that  the  shroud  tube  itself  was  below  bath 
temperature.  The  net  effect  of  the  slug  formation  was  to  cause  the 
combustor  ullage  region  pressure  to  exceed  that  of  the  shroud  tube 
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Table  2 


Summary  of  Air-Cooled  Combustor  Tests 


Run 

Number 

Run 

Length 

(hours) 

Oxidizer  Flow 
(kg/hr) 

Coolant  Flow 
(kg/hr) 

Maximum 

Bath 

Temperature 

(°K) 

Maximum 

Power 

Level 

(kW) 

PPC-1 

2.4 

2.26 

22.0-84.9 

1208 

10.8 

PPC-2 

2.6 

0.90-2.27 

31.8-43.0 

1213 

10.8 

PPC-3 

1.4 

2.27 

20.1-48.8 

1225 

10.8 

PPC-4 

3.8 

0.75-2.25 

0.0 

1210 

10.7 

PPC-5 

2.5 

0.75-1.88 

16.4-37.7 

1172 

8.9 

PPC-6 

0.9 

0.76 

0.0 

1161 

3.6 

PPC-7 

2.8 

0.56-1.13 

0.0 

1238 

5.4 

PPC-8 

5.2 

0.47-1.13 

0.0 

1255 

5.4 

PPC-9 

6.3 

1.31-4.96 

18.0-109. 

1205 

23.6 

PPC-10 

10.0 

1.11-5.29 

11.0-127. 

1205 

25.2 
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during  the  period  of  time  the  slug  was  moving  through  the  shroud 
tube.  The  pressure  mismatch  upset  the  hydrostatic  balance  of  the 
combustor  trap  which  resulted  in  the  loss  of  the  product  liquid  in 
the  trap  and  ultimate  trap  failure.  Normal  trap  operation  had 
occurred  in  both  tests  prior  to  the  slug  formation. 

Several  factors  existed  which  made  the  air-cooled  combustor 
more  sensitive  to  slug  formation  than  the  radiative  combustor,  which 
never  showed  signs  of  this  problem.  The  factors  include:  a)  a 
greater  argon  flow  into  the  combustor  ullage  region;  b)  a  longer 
shroud  tube  length;  and  c)  a  smaller  product  liquid  volume.  The 
greater  argon  flow  rate  resulted  in  a  larger  rate  of  product  outflow 
during  the  pressurization  period,  and  the  longer  shroud  tube 
increased  the  period  as  the  slug  had  to  move  a  greater  distance. 

The  two  effects  were  magnified  for  the  air-cooled  combustor,  due  to 
the  smaller  product  volume  which  could  be  outflowed  before  the  trap 
failed. 

Tests  PPC-3  and  PPC-4  were  run  with  a  spot  heater  installed 
in  the  crossover  port  area  to  maintain  the  temperature  of  the  region 
close  to  bath  temperature.  Both  tests  were  terminated  when  failures 
of  the  heater  occurred.  It  was  found  that  the  spot  heater  could  not 
maintain  the  top  of  the  combustor  at  bath  temperature  when  coolant 
air  was  applied  to  the  system,  and  the  condition  led  to  failure  of 
the  trap  as  described  above  during  test  PPC-5.  The  initial  design  of 
the  combustor  did  not  permit  the  proper  installation  of  a  spot  heater. 
Alternate  methods  of  installation  were  tried  on  tests  PPC-6,  PPC-7 
and  PPC-8,  all  of  which  forced  the  tests  to  be  terminated  when  the 
heaters  failed.  The  three  tests  were  run  at  low  power  levels 
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which  permitted  the  argon  flow  to  the  system  through  the  coaxial 
injector  to  be  reduced,  and  eliminated  the  need  for  coolant  air  flow. 
No  problems  were  encountered  with  the  trap  during  the  tests,  but  it 
was  evident  that  a  design  change  would  be  needed  before  the  system 
could  be  operated  at  higher  power  levels. 

Bath  liquid  samples  were  obtained  during  tests  PPC-2  through 
PPC-8.  The  samples  were  analyzed  following  each  test,  and  excellent 
agreement  was  obtained  between  samples  obtained  from  the  different 
tests . 

The  combustor  was  modified  prior  to  the  final  two  tests  to 
the  configuration  described  in  Section  4.2.  Tests  PPC-9  and  PPC-10 
were  conducted  to  complete  the  remaining  test  objective  of  measuring 
the  thermal  performance  of  the  system. 

Power  levels  of  6.2,  12.8,  19.4,  21.3,  and  23.4  kW  were 
examined  during  test  PPC-9  by  operating  at  each  condition  for  one 
hour.  The  oxidizer  flow,  coolant  air  flow,  bath  temperature,  exhaust 
air  temperature,  and  mean  cycle  load  for  the  test  are  plotted  in 
Figure  13  as  a  function  of  time.  The  coolant  air  flow  was  adjusted 
at  each  test  condition  to  maintain  the  bath  temperature  at  1197+3°K. 

A  total  of  19.14  kg  of  oxidizer  and  7.25  kg  of  fuel  were  consumed 
during  the  test.  The  test  was  terminated  when  all  test  objectives 
were  met. 

Additional  performance  data  were  obtained  during  test  PPC-10 


by  operating  for  one  hour  periods  at  power  levels  of  7.2,  8.2,  10.1, 
12.5,  13.9,  15.3,  20,5,  and  23.9  kW.  The  test  was  unique  in  that 
both  the  longest  run  time,  10  hours,  and  highest  power  output,  25.2 
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kW,  were  achieved  during  this  test.  Figure  14  presents  the  oxidizer 
flow,  coolant  air  flow,  bath  temperature,  exhaust  air  temperature,  and 
mean  cycle  load  for  the  test  plotted  versus  time.  The  coolant  air 
flow  was  adjusted  at  each  condition  to  maintain  the  same  bath 
temperature  examined  in  test  PPC-9 ,  119 7+3 °K.  A  total  of  26.4  kg 
of  oxidizer  and  10,0  kg  of  fuel  were  consumed  during  the  test,  and 
the  test  was  terminated  when  all  the  objectives  were  met.  The 
performance  data  obtained  during  test  PPC-10  agreed  well  with  the 
data  obtained  during  test  PPC-9. 

5.1.3  Nature  of  the  Reaction  Products  The  reaction 
products  produced  during  the  radiative  and  air-cooled  combustor 
tests  were  collected  by  a  dry-tank  collection  method  where  the 
products  are  exhausted  into  a  tank  originally  empty  at  the  start  of 
the  test.  Figure  15  is  a  photograph  of  the  product  material  collected 
during  test  PPC-9.  The  products  formed  a  stalagmite  formation 
extending  from  the  bottom  of  the  tank  to  a  height  of  approximately 
one  meter.  The  mass  of  the  product  shown  in  Figure  15  was 
approximately  25  kg.  The  product  formation  produced  during  test 
PPC-9  was  typical  of  the  formations  produced  during  other  tests. 

An  alternate  method  of  product  collection  is  the  wet-tank 
collection  method  where  the  products  are  exhausted  into  the  fuel 
tank  containing  the  molten  fuel.  Such  a  method  is  possible  because 
the  volume  of  the  reaction  products  is  less  than  the  volume  of  the 
fuel.  The  method  would  result  in  a  reduction  in  overall  system 
volume.  The  wet-tank  collection  method  was  not  investigated  in  an 
actual  combustor,  but  was  investigated  using  the  apparatus  described 
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Figure  14  Air-Cooled  Combustor  Run  Data  (Test  PPC-10) 
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Photograph  of  Reaction  Products  Figure  16  Photograph  of  Wet-Tank  Reaction 

Products 
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in  Appendix  F.  The  form  of  the  product  collected  with  a  wet-tank 
collection  method  is  shown  as  a  photograph  in  Figure  16.  The  products 
formed  spherical  balls  approximately  0.9  cm  in  diameter.  The  packing 
efficiency  of  the  balls  was  approximately  50%  on  a  volume  basis. 

5.2  Liquid  Phase  Solubility  Results 

The  theoretical  and  experimental  solubility  results  for  the 
fuel-rich  phase  are  shown  in  Figure  17  as  mole  fraction  lithium 
fluoride  and  lithium  sulfide  plotted  versus  temperature  with  the 
remainder  of  the  phase  consisting  of  pure  lithium.  The  experimental 
data  is  presented  in  Table  3.  The  data  taken  from  different  combustor 
runs  are  distinguished  by  different  symbols  in  the  figure,  and  it  is 
seen  that  the  agreement  between  the  different  data  is  good.  The 
theoretical  results  obtained  from  the  thermodynamic  model  described 
in  Chapters  II  and  III  are  shown  as  solid  lines  in  Figure  17  for  the 
case  employing  the  high  value  for  the  heat  of  fusion  of  lithium  sulfide, 
49.4  kJ/mole,  and  as  a  dashed  line  for  the  case  employing  the  low 
value  for  the  heat  of  fusion,  7.24  kJ/mole.  The  results  for  the  two 
cases  agreed  within  0.6%  for  the  lithium  fluoride  component.  The 
agreement  between  the  theoretical  predictions  and  the  experimental 
data  is  good.  The  maximum  error  in  the  lithium  fluoride  solubility 
prediction  is  approximately  25%  and  occurs  at  the  low  temperature  end 
of  the  temperature  range  examined.  The  maximum  error  in  the  lithium 
sulfide  solubility  prediction  is  approximately  39%  for  the  high  heat 
of  fusion  case,  and  occurs  at  the  high  temperature  end  of  the 
temperature  range.  The  error  for  the  low  heat  of  fusion  case 
remains  relatively  constant  at  a  value  of  approximately  50%.  Both 
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Table  3 


; 

i 

|  Fuel-Rich  Liquid  Solubility  Data 


Test 

Number 

Temperature 

<°K) 

Lithium  Mole 
Fraction 

Lithium  Fluoride 
Mole  Fraction 

Lithium  Sulfide 
Mole  Fraction 

PPC-4 

1116 

0.980 

0.0172 

0.00233 

PPC-8 

1116 

0.977 

0.0196 

0.00309 

PPC-8 

1136 

0.976 

0.0209 

0.00334 

PPC-7 

1145 

0.974 

0.0228 

0.00357 

PPC-5 

1147 

0.971 

0.0255 

0.00364 

PPC-8 

1155 

0.973 

0.0236 

0.00372 

PPC-8 

1161 

0.970 

0.0260 

0.00432 

PPC-7 

1166 

0.970 

0.0255 

0.00401 

PPC-8 

1183 

0.961 

0.0338 

0.00553 

PPC-4 

1186 

0.974 

0.0222 

0.00336 

PPC-4 

1200 

0.962 

0.0333 

0.00502 

PPC-4 

1202 

0.970 

0.0259 

0.00423 

PPC-7 

1214 

0.959 

0.0352 

0.00558 

PPC-8 

1214 

0.958 

0.0359 

0.00554 

PPC-7 

1238 

0.948 

0.0448 

0.00717 

PPC-7 

1238 

0.944 

0.0490 

0.00714 

PPC-8 


1244 


0.958 


0.0362 


0.00630 
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the  experimental  data  and  the  theoretical  predictions  show  a  gradual 
increase  in  lithium  fluoride  and  lithium  sulfide  solubilities  in  the 
fuel-rich  phase  as  the  temperature  is  increased. 

Figure  18  compares  the  theoretical  and  experimental  solubility 
results  for  the  product-rich  phase,  and  Table  4  presents  a  summary  of 
the  experimental  data.  The  lithium  solubility  is  plotted  as  mole 
fraction  versus  temperature.  The  ratio  of  mole  fraction  lithium 
sulfide  to  mole  fraction  lithium  fluoride  is  plotted  versus  temperature 
since  the  solubility  ratio  of  these  two  components  is  controlled  by 
the  reaction  stoichiometry.  The  agreement  between  the  experimental 
data  taken  from  different  combustor  runs,  distinguished  by  different 
symbols  in  the  figure,  is  seen  to  be  good.  The  theoretical  results 
computed  from  the  two  different  heat  of  fusion  cases  are  identical 
within  the  accuracy  of  the  plot,  and  are  shown  as  solid  lines. 

The  lithium  solubility  computed  from  the  low  heat  of  fusion  case 
was  approximately  4%  lower  than  the  high  heat  of  fusion  case  at 
1110°K  and  2%  higher  at  1250° K.  The  theoretical  predictions  are 
extremely  good  for  the  product-rich  phase.  The  lithium  solubility 
is  predicted  to  an  accuracy  of  5-8%  over  the  temperature  range  examined. 
The  lithium  fluoride-lithium  sulfide  mole  fraction  ratio  is 
theoretically  exact  and  any  disagreement  is  caused  by  scatter  of  the 
experimental  data.  The  lithium  solubility  in  the  product-rich  phase 
was  found  to  be  very  low,  less  than  0.9%  on  a  mole  fraction  basis, 
in  the  temperature  range  examined,  and  was  found  to  increase  with 
temperature  as  expected. 

The  scatter  in  the  data  presented  in  Figures  17  and  18  could 
be  caused  by  the  following:  a)  nonequilibrium  of  the  phases  when 
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Table  4 


Product-Rich  Liquid  Solubility  Data 


Test 

Number 

Temperature 

(°K) 

Lithium  Mole 
Fraction 

Lithium  Fluoride 
Mole  Fraction 

Lithium  Sulfide 
Mole  Fraction 

PPC-8 

1119 

0.00200 

0.856 

0.141 

PPC-8 

1130 

0.00274 

0.864 

0.133 

PPC-7 

1144 

0.00425 

0.865 

0.131 

PPC-5 

1144 

0.00235 

0.851 

0.147 

PPC-7 

1147 

0.00444 

0.859 

0.136 

PPC-8 

1158 

0.00274 

0.864 

0.133 

PPC-8 

1166 

0.00278 

0.857 

0.140 

PPC-8 

1172 

0.00346 

0.858 

0.138 

PPC-7 

1172 

0.00381 

0.857 

0.139 

PPC-5 

1180 

0.00380 

0.854 

0.142 

PPC-7 

1200 

0.00448 

0.853 

0.142 

PPC-3 

1205 

0.00637 

0.868 

0.125 

PPC-2 

1205 

0.00675 

0.849 

0.144 

PPC-8 

1227 

0.00662 

0.848 

0.145 

PPC-8 

1230 

0.00638 

0.850 

0.145 

PPC-8 

1233 

0.00665 

0.852 

0.141 

the  bath  samples  were  taken;  b)  the  sample  composition  not 
representative  of  the  bulk  liquid-composition;  c)  errors  in  temperature 
measurement;  and  d)  errors  in  determining  the  bath  sample  composition. 
As  described  in  Section  4.5.3,  precautions  were  taken  during 
the  operation  of  the  combustor  to  ensure  that  equilibrium  conditions 
were  established  in  the  combustor.  A  quantitative  estimate  of 
nonequilibrium  effects  is  not  available,  but  it  is  noted  that  the 
rate  of  phase  saturation  is  enhanced  by  the  very  low  solubilities 
of  the  minor  components  involved  in  each  phase.  The  accuracy  of 
the  thermocouples  is  listed  by  the  manufacturer  to  be  within 
approximately  +9°K  in  the  temperature  range  of  interest.  The 
measurement  error  introduced  should  be  consistent  from  test  to  test 
because  the  thermocouples  used  in  the  different  tests  were  taken 
from  the  same  production  lot.  The  accuracy  with  which  the  recorded 
temperatures  could  be  read  from  the  chart  was  estimated  to  be  +2°K. 

It  was  concluded  that  the  measurement  accuracy  of  the  thermocouples 
was  well  within  the  manufacturer's  tolerance  since  the  product 
eutectic  temperature,  1065°K,  and  the  fuel  solidification  temperature, 
453°K,  were  consistently  measured  to  within  the  accuracy  of  the 
recording  chart  by  observing  thermal  arrests  during  post-test  cool¬ 
down  cycles.  Errors  introduced  by  the  wet-test  chemical  analysis 
procedures  were  found  to  be  small  by  testing  the  procedures  against 
prepared  samples  of  known  composition.  The  lithium  fluoride  and 
lithium  sulfide  concentrations  in  the  fuel-rich  and  product-rich 
phases  were  found  to  be  within  an  accuracy  of  +3%.  The  lithium 
concentration  in  the  product-rich  phase  was  determined  to  within 
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an  accuracy  of  +8%.  The  greatest  source  of  error  could  not  be 
evaluated  quantitatively,  but  was  believed  to  be  caused  by  the  bath 
samples  not  representing  the  bulk  phase  composition.  The  composition 
change  could  occur  after  sample  removal  by  small  quantities  of 
lithium  being  lost  through  the  ball  valve/sample  cup  seal  in  the 
case  of  the  fuel-rich  samples ,  or  by  contamination  of  the  samples 
with  extraneous  lithium  in  the  case  of  the  product-rich  samples. 

5 . 3  Liquid  Phase  Density  Results 

The  theoretical  and  experimental  densities  for  both  the  product- 
rich  and  fuel-rich  liquid  phases  are  compared  in  Figure  19.  The 
experimental  data  are  presented  in  Table  5.  The  theoretical  results 
for  the  two  cases  examined  in  the  analysis  (involving  different  heats 
of  fusion  for  the  lithium  sulfide  component)  are  the  same  within  the 
accuracy  of  the  plot.  The  largest  difference  between  the  two  cases 
is  approximately  0.7%  for  the  fuel-rich  phase  density  at  1110°K. 

The  density  of  the  fuel-rich  phase  was  found  experimentally  to  be 
relatively  constant  in  the  temperature  range  examined.  The  theoretical 
result,  shown  as  the  solid  line  in  Figure  19,  predicts  a  slight 
increase  in  density  with  temperature  and  gives  values  approximately 
6%  lower  than  the  experimental  data  in  the  low  temperature  region. 

The  agreement  in  the  higher  temperature  region  is  within  approximately 
2%.  The  lack  of  temperature  dependence  of  fuel-rich  liquid  phase 
density  is  caused  by  a  balance  between  the  decrease  in  the  density  of 
pure  lithium  and  the  increased  solubility  of  the  more  dense  product 
species  with  increasing  temperature.  The  product-rich  liquid  density 
was  found  to  decrease  with  increasing  temperature.  The  data  agree 
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Figure  19  Bath  Densities 


Table  5 


Bath  Liquid  Density  Data 


Temperature 

(°K) 


Fuel-Rich  Liquid  Density 
_ (kg/m3) _ 


Product-Rich  Liquid  Density 
_ (kg/m3) _ 
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with  the  theoretical  prediction,  shown  as  the  solid  line  in  Figure  19, 
to  within  +2%  except  at  the  highest  temperature  where  the  error 
reaches  a  maximum  of  approximately  5%. 

5.4  Gas  Phase  Partial  Pressure  Results 

The  theoretical  gas  phase  partial  pressures  are  presented  in 
Figure  20  as  the  sum  of  the  partial  pressures  of  various  dimer  and 
polymeric  species  plotted  versus  temperature.  The  predictions  are 
effectively  unchanged  by  the  value  of  the  heat  of  fusion  of  lithium 
sulfide  used  in  the  model;  the  predictions  for  the  lithium  and  lithium 
fluoride  species  agreeing  within  0.1%  and  the  predictions  for  the 
lithium  sulfide  species  agreeing  within  0.25%.  The  partial 
pressures  of  the  individual  species  were  not  measured  experimentally, 
but  the  total  pressure  was  measured  during  the  course  of  a  related 
study  [46],  The  method  used  was  to  measure  the  oxidizer  inlet 
pressure  to  an  unvented  wick-type  pot  combustor.  The  measured 
pressure  was  equal  to  the  mixture  vapor  pressure  because  the  injector 
was  unchoked,  and  the  combustion  was  occurring  off  the  surface  of 
the  wick  as  a  diffusion  flame.  The  experimental  total  pressure 
results  are  compared  to  the  theoretical  predictions  in  Figure  21. 

The  largest  differences  of  10-15%  occur  in  the  low  temperature  region 
where  the  experimental  uncertainties  were  the  greatest.  The  results 
indicate  that  the  reaction  chamber  pressure  in  unvented  combustors 
can  vary  approximately  an  order  of  magnitude  in  the  typical  operating 
temperature  range,  but  in  all  cases  is  very  low  (0.035  bar  at  1250°K). 
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5.5  Oxidizer  Utilization  Efficiency  Results 

Oxidizer  utilization  efficiency  tests  were  conducted 
following  the  procedure  described  in  Appendix  E.  The  injector  was 
submerged  at  depths  of  5.0  to  7.6  cm.  and  had  oxidizer  and  shroud 
argon  flows  of  approximately  2.5  and  0.63  kg/hr,  respectively.  The 
tests  were  conducted  in  a  vented  pot-type  combustor.  The  results 
indicated  the  presence  of  small  amounts  of  unreacted  oxidizer  in  the 
argon  flow  leaving  the  combustor  during  the  initial  start-up  transient 
when  the  bath  temperature  was  in  the  range  of  810  to  870°K.  No 
unreacted  oxidizer  was  detected  during  subsequent  operation.  The 
oxidizer  utilization  efficiency  during  the  start-up  transient  was 
greater  than  99.5%.  The  results  confirmed  that  the  oxidizer 
utilization  efficiency  is  100%  during  normal  operating  conditions 
of  the  experimental  apparatus. 

5.6  Thermal  Performance  Results 

5.6.1  Energy  Release.  The  thermal  performance  data  obtained 
during  tests  PPC-9  and  PPC-10  are  plotted  in  Figure  22  as  heat 
exchanger  load  versus  oxidizer  flow.  The  average  net  parasitic  heat 
loss  from  the  apparatus  is  given  by  the  negative  intercept  on  the 
heat  exchanger  load  axis,  and  includes  the  sum  of  all  losses  from 
the  system  and  the  electrical  input  to  the  system  through  the  sump 
guard  heaters.  The  value  of  2.98  kW  was  determined  by  a  least 
squares  linear  fit  to  the  data  and  therefore  is  average  value  for 
the  entire  group  of  data.  The  slope  of  the  least  squares  fit  yields 
an  energy  release  rate  of  4.99  kW-hr/kg  of  oxidizer  at  a  bath 
temperature  of  1197+3°K,  a  fuel  inlet  temperature  of  533°K,  and  an 


95 


oxidizer  temperature  of  305°K.  The  theoretical  values  for  these 
conditions  are  4.71  kW-hr/kg  and  4.79  kW-hr/kg  for  the  high  and 
low  lithium  sulfide  heats  of  fusion  cases,  respectively.  The 
experimental  value  is  5.9%  higher  than  the  high  heat  of  fusion  value 
and  4.2%  higher  than  the  low  heat  of  fusion  value.  The  small  error 
largely  results  from  the  assumption  that  the  parasitic  heat  loss 
remained  constant  with  load.  In  fact,  the  parasitic  heat  loss 
decreased  with  increased  load.  The  value  of  the  parasitic  heat 
loss  at  each  load  condition  could  not  be  determined,  and  the  use  of 
an  average  value  caused  the  slope  of  the  curve  to  be  incorrectly  too 
high. 

The  difference  between  the  theoretical  energy  release  rates 
is  caused  by  the  enthalpy  of  liquid  lithium  sulfide  being  approximately 
10%  higher  for  the  high  heat  of  fusion  case  than  for  the  low  heat  of 
fusion  case.  It  is  impossible  to  determine  from  the  thermal 
performance  data  which  is  the  correct  value  for  the  heat  of  fusion 
of  lithium  sulfide,  as  errors  introduced  by  the  estimated  heat 
capacity  data  are  expected  to  be  approximately  of  the  same  order  of 
magnitude  as  the  changes  caused  by  the  differences  in  the  heats 
of  fusion.  Additionally,  the  uncertainty  in  the  experimental  value 
dominates  the  small  difference  between  the  two  theoretical  values. 

The  contribution  to  the  product-rich  liquid  enthalpy  from  the  excess 
enthalpy  of  mixing  is  on  the  order  of  0.1%  and  therefore  effects 
from  nonideal  mixing  are  small  from  the  energy  release  standpoint. 

The  recommended  value  for  the  energy  release  rate  at  the  given 
operating  conditions  is  the  average  of  the  two  results  or  4.75  kW-hr/kg 


of  oxidizer. 
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5.6.2  Energy  Utilization.  The  energy  utilization  efficiency 
of  the  air-cooled  combustor  is  considered  in  Figure  23.  The  figure 
presents  an  energy  flow  diagram  for  the  combustor  under  conditions 
where  the  power  to  the  load  is  25  kW  and  the  bath  temperature,  fuel 
inlet  temperature  and  oxidizer  temperature  are  1197,  533,  and  305°K, 
respectively.  The  losses  from the  system  associated  with  parasitic 
heat  loss  and  product  enthalpy  were  found  to  equal  3.7  and  5.57 
kW,  respectively.  The  parasitic  heat  loss  was  determined  by 
extrapolating  the  data  at  the  22.5  kW  load  condition  and  using  the 
average  theoretical  energy  release  rate  of  4.75  kW-hr/kg  to  obtain 
the  net  parasitic  heat  loss.  The  power  input  from  the  electrical 
guard  heaters  of  1.9  kW  was  estimated  from  knowledge  of  the  heater 
power  outputs  and  the  controller  cycle  periods.  The  electrical  power 
was  added  to  the  net  parasitic  heat  loss  to  obtain  the  total 
parasitic  heat  loss  from  the  system.  It  is  noted  that  the  net 
parasitic  heat  loss  from  the  system  at  the  extrapolated  25  kW 
power  level  was  approximately  40%  lower  than  the  average  value  for  all 
actual  load  conditions.  The  product  enthalpy,  electrical  energy 
input  to  the  fuel,  and  the  standard  heat  of  reaction  were  computed  from 
the  theoretical  model  using  the  average  energy  release  rate. 

In  an  actual  submersible  application  the  product  enthalpy 
would  not  be  lost  from  the  system,  but  could  be  used  as  a  source  for 
vehicle  heating  and  as  energy  needed  to  melt  the  fuel.  Additionally, 
the  electrical  guard  heaters  could  be  eliminated  by  designing  the 
sump  regions  with  additional  emphasis  on  the  thermal  management 
problems  associated  with  these  regions. 
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As  illustrated  in  Figure  23,  approximately  79.2%  of  the  standard 
heat  of  reaction  at  298°K  was  transferred  to  the  load  at  a  bath 
temperature  of  1197°K.  The  effects  of  bath  temperature  on  system 
performance  are  discussed  in  the  following  section. 

5.6.3  Bath  Temperature  Effects  on  System  Performance.  Figure 
24  presents  fuel  utilization  efficiency,  thermal  efficiency,  bath 
liquid  density  ratio,  and  total  mixture  pressure  as  a  function  of 
temperature.  The  results  are  the  average  of  the  values  computed 
from  the  two  cases  employing  different  values  for  the  heat  of  fusion 
of  lithium  sulfide  examined  in  the  theoretical  model.  The  thermal 
efficiency  curve  was  computed  assuming  no  parasitic  heat  loss  from 
the  system  and  with  fuel  and  oxidizer  inlet  temperatures  of  533  and 
305 °K,  respectively.  Additionally,  the  effects  of  fuel  carryover 
in  the  argon  exhaust  from  the  reaction  chamber  and  the  energy  loss 
resulting  from  the  flow  were  not  included. 

The  bath  liquid  density  ratio  given  by  the  ratio  of  the 
density  of  the  product-rich  liquid  to  the  density  of  the  fuel-rich 
liquid  decreases  with  increasing  bath  temperature  from  a 
value  of  3.52  at  1110°K  to  a  value  of  1.00  at  the  computed  consolute 
temperature  of  1638°K.  At  temperatures  above  the  consolute  temperature 
the  product-rich  and  fuel-rich  phases  are  miscible  and  removal  of  the 
reaction  products  from  the  reaction  chamber  is  not  possible.  The 
decrease  in  the  liquid  density  ratio  with  increasing  temperature 
indicates  that  the  combustor  trap  geometry  needed  to  maintain  the 
fuel-rich  and  product-rich  liquid  levels  at  specified  values  is 
temperature  dependent.  For  the  trap  design  employed  in  the  combustor 
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used  in  the  present  study,  the  inner  trap  tube  would  have  to  be 
lengthened  for  operation  at  a  higher  bath  temperature. 

The  total  vapor  pressure  of  the  mixture  was  found  to  increase 
significantly  with  increasing  bath  temperature.  The  total  pressure 
increases  from  0.00565  bar  at  1110°K  to  0.98  bar  at  1638°K.  The 
vapor  pressure  of  the  mixture  becomes  the  internal  pressure  of  the 
reaction  chamber  when  an  unvented  configuration  is  used.  In  this 
case,  the  combustor  trap  system  must  be  designed  with  the  pressure 
forces  taken  into  consideration.  For  a  vented  combustor  such  as  the 
one  used  in  the  present  study,  the  increased  vapor  pressure  represents 
an  increased  mass  flow  of  fuel  and  product  species  from  the  combustor 
for  a  given  argon  flow. 

The  thermal  efficiency  is  defined  as  the  ratio  of  the  energy 
released  to  the  load  to  the  standard  heat  of  reaction  at  298°K. 

The  thermal  efficiency  decreases  from  a  value  of  86.4%  at  a  bath 
temperature  of  1110° K  to  a  value  of  67.0%  at  a  bath  temperature  of 
1638°K.  The  decrease  is  caused  by  the  products  being  exhausted  at  a 
higher  temperature  and  therefore  with  a  higher  enthalpy. 

The  fuel  utilization  efficiency  is  defined  as  the  ratio  of  the 
rate  of  fuel  consumption  by  the  chemical  reaction  to  the  total  fuel 
flow  rate  to  the  reaction  chamber.  Fuel  utilization  efficiencies 
of  99.8%  and  73.0%  are  achieved  at  bath  temperatures  of  1110  and 
1638°K,  respectively.  The  decrease  is  caused  by  the  increased 
solubility  of  the  fuel  in  the  products  exhausted  from  the  reaction 
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CHAPTER  VI 


SUMMARY  AND  CONCLUSIONS 


2. 


6.1  Summary 

The  overall  objective  of  the  present  study  was  to  investigate 
analytically  and  experimentally  the  characteristics  of  a  steadily 
operating  liquid  metal  thermal  energy  source  employing  the  lithium- 
sulfur  hexafluoride  reactant  combination.  The  specific  objectives 
of  the  study  were  as  follows: 

1.  Experimentally  measure  bath  liquid  densities,  bath  liquid 
compositions,  and  thermal  performance  for  temperatures 
typical  of  metal  combustor  operation. 

Develop  a  thermodynamic  model  for  the  bath  liquids  and 
compare  the  model  predictions  to  the  measurements. 

Examine  continuous  operation  of  a  complete  steady 
combustor  system  with  varying  thermal  load  requirements 
for  extended  periods  of  time. 

A  steady  metal  combustor  apparatus  which  permitted  the 
necessary  measurements  to  be  made  was  designed  and  fabricated. 
Steady-state  conditions  were  established  by  continuously  supplying 
fuel  and  oxidizer  to  the  reaction  chamber  while  continuously  removing 
the  reaction  products  through  a  trap  arrangement  in  the  bottom  of  the 
chamber.  The  oxidizer  was  injected  into  the  molten  bath  through  a 
coaxial  injector,  and  fuel  was  supplied  with  an  argon  lift  pump. 

The  apparatus  was  operated  in  the  temperature  range  of  1110  to  1255°K 
with  power  outputs  ranging  from  0.0  to  25.2  kW. 
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Procedures  were  developed  to  obtain  bath  liquid  samples  from 
the  reaction  chamber,  and  to  analyze  the  samples  using  wet-test 
chemical  analysis  methods.  Procedures  were  also  developed  for 
obtaining  bath  liquid  density  data.  Bath  liquid  solubility  data  were 
obtained  in  the  temperature  range  of  1116  to  1244 °K,  and  bath  liquid 
density  data  were  obtained  in  the  temperature  range  of  1144  to  1280°K, 

Thermal  performance  data  were  generated  at  a  bath  temperature 
of  1197+3°K  over  a  range  of  load  conditions.  The  data  yielded  an 
energy  release  rate  of  4.99  kW-hr/kg  of  oxidizer  which  is  in  good 
agreement  with  the  theoretical  value  obtained  from  the  thermodynamic 
model  of  4.75  kW-hr/kg. 

System  operation  was  demonstrated  for  extended  periods  with  the 
longest  test  being  ten  hours  in  duration.  The  capability  of  the  system 
to  operate  under  variable  load  range,  idling  and  restart  conditions 
was  also  demonstrated.  No  corrosion  problems  with  the  type  316 
stainless  steel  components  were  observed. 

A  thermodynamic  model  for  the  reaction  chamber  contents  was 
developed  using  the  van  Laar  expression  for  the  excess  Gibbs  free 
energy  of  mixing  in  the  liquid  phase  and  an  ideal  solution 
assumption  for  the  gas  phase.  The  model  permitted  the  calculation 
of  liquid  phase  solubilities,  liquid  phase  densities,  gas  phase 
partial  pressures,  and  thermal  performance. 

Experimentally  determined  bath  liquid  solubility  data, 
bath  liquid  density  data,  and  thermal  performance  were  compared 
to  the  theoretical  predictions.  It  was  found  that  the  model 
adequately  predicted  the  thermodynamic  characteristics  of  the  system. 
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6.2  Conclusions 

The  major  conclusions  of  the  study  are  as  follows: 

1.  The  thermodynamic  model  for  the  ternary  immiscible  liquid 
mixture  composed  of  lithium,  lithium  fluoride,  and 
lithium  sulfide  examined  in  this  study  predicted  the 
solubilities  of  the  various  components  within  an 
accuracy  of  5-50%  indicating  the  van  Laar  expression  for 
the  excess  Gibbs  free  energy  of  mixing  used  in  the  model 
adequately  represents  the  mixing  characteristics  of  the 
system.  Additionally,  the  experimentally  measured 
liquid  densities,  mixture  vapor  pressures  and  thermal 
performance  agreed  within  5%  of  the  thermodynamic  model 
predictions. 

2.  A  steadily  operating  metal  combustor  was  operated  for 
periods  on  the  order  of  ten  hours  under  conditions  of 
continuous  fuel  addition  and  product  outflow,  variable 
load  range,  idling  and  restart.  These  operating  periods 
were  not  limited  by  the  combustor,  and  longer  periods  of 
operation  appear  to  be  feasible. 

3.  The  thermal  efficiency  of  the  system  decreases  with 
temperature  from  a  value  of  86.4%  at  a  bath  temperature 
of  1110°K  to  a  value  of  84.0%  at  a  bath  temperature  of 
1250°K. 

4.  Fuel  utilization  efficiency  decreases  with  increasing  bath 
temperature  from  a  value  of  99,8%  at  a  bath  temperature  of 
1110°K  to  a  value  of  99.2%  at  a  bath  temperature  of  1250°K. 
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5.  Visual  inspection  of  the  type  316  stainless  steel 

components  used  in  the  experimental  apparatus  revealed 
no  corrosion  problems,  indicating  that  this  material  is 
suitable  for  combustor  fabrication.  Sample  contamination 
was  found  during  the  chemical  analysis  of  fuel-rich  liquid 
samples  extracted  with  type  316  stainless  steel  sample 
cups  which  was  not  found  in  samples  extracted  with  nickel 
200  alloy  cups.  The  contamination  was  believed  to  be 
caused  by  the  products  of  surface  corrosion  of  the 
stainless  steel  by  lithium  sulfide  since  similar 
contamination  was  found  during  lithium-lithium 
sulfide  solubility  tests  at  temperatures  above  1100°K, 
These  findings  suggest  low  levels  of  corrosion  that  are 
not  apparent  during  visual  inspection  of  components 
following  test  periods  in  the  range  of  the  current 
investigation,  but  which  could  be  a  factor  during  longer 
periods  of  operation. 

6.3  Recommendations  for  Further  Study 

It  is  desirable  to  have  additional  solubility  data  for  the 
lithium-lithium  sulfide  and  lithium  fluoride-lithium  sulfide  binary 
systems.  Additional  data  over  a  broader  temperature-concentration 
range  would  permit  more  confident  evaluation  of  the  empirical 
parameters  associated  with  these  binary  pairs.  Of  particular 
importance  is  data  in  the  lithium  sulfide-rich  region  of  either  phase 
diagram  since  such  data  could  be  used  to  infer  the  heat  of  fusion  of 


lithium  sulfide. 
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The  thermodynamic  model  developed  in  this  study  can  be 
applied  to  the  analysis  of  other  combustor  concepts.  Several  new 
combustor  designs  have  been  proposed  which  need  to  be  examined  both 
analytically  and  experimentally  [46]. 
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APPENDIX  A 


COMPUTER  PROGRAM  AND  NEWTON-RAPHSON  METHOD  USED 
TO  SOLVE  THE  THERMODYNAMIC  MODEL  EQUATIONS 

The  set  of  thermodynamic  model  equations  are  not  easily  solved 
because  of  the  nonlinear  form  of  Equations  (2.43)  after  substitution 
of  Equations  (2.19)  for  the  component  activities.  The  resulting 
equations  were  simplified  somewhat  by  taking  the  natural  logarithm  of 
both  sides. 

The  Newton-Raphson  method  was  used  to  obtain  the  solution 
to  the  set  of  nonlinear  equations.  A  change  of  variables  was 
necessitated  because  the  existence  of  the  natural  logarithm  of  the 
mole  fractions  resulted  in  procedure  failure  if  the  computer 
iterated  to  a  negative  mole  fraction  during  the  course  of  the  search. 
To  eliminate  the  problem  logarithmic  variables  were  employed  as 
follows: 


“  In  x^ 

(A.l) 

-  In  y± 

(A.  2) 

where  x^  and  y^  are  the  mole  fractions  of  the  various  components  in  the 
two  immiscible  phases. 

In  terms  of  the  variables  defined  by  Equations  (A.l)  and  (A. 2), 
Equations  (2.43),  (2.44),  (2.45),  and  (2.49)  become 
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Equations  (A. 3)  to  (A. 8)  were  linearized  by  choosing  estimates 
of  the  variables,  z^Q  and  w^q,  and  writing 


z.mz.  +  z. ' 
i  io  i 


w  ■  w.  +  W  ' 
i  io  i 


(A. 9) 
(A. 10) 


where  z^'  and  w^'  are  small  correction  terms  to  the  original  estimate. 

The  linearized  equations  can  be  obtained  either  by  direct 
substitution  of  Equations  (A. 9)  and  (A. 10)  into  Equations  (A. 3)  to 
(A. 8)  and  simplification,  or  by  expanding  each  equation  in  a  Taylor 
series  about  the  initial  estimates.  Using  the  Taylor  series  method, 
the  equations  can  be  expressed  as  functions  of  the  six  variables 
as  follows 


f(zi,  w±)  -  C 


(A. 11) 


where  C  represents  a  constant  characteristic  of  a  given  equation.  If 
higher  order  terms  are  neglected,  the  equations  can  be  expressed  as 

3 


!(*i*  wi)  “  f(zio*  Wio)  +  J1*i’  (zi0’  Wio) 


+  5lWi'  "i°)  “  0 


(A. 12) 


The  resulting  expressions  are  arranged  in  matrix  form  as  follows 


tv 


w 


L  i  J 


-  G, 


(A. 13) 


The  computer  program  used  to  solve  Equations  (A. 13)  is  presented 
below  with  the  matrix  coefficients  and  being  denoted  by  A(i,j) 
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and  R(i)  in  the  program.  The  equations  were  solved  for  z^'  and  w^ 
for  a  given  iteration  using  the  Gaussian  elimination  method.  The 
original  estimates  were  updated  following  each  iteration  by  adding 
z^'  and  w^'  to  ziQ  and  w^Q,  respectively,  and  the  method  was 
repeated  until  all  z '  and  w^'  were  reduced  to  within  desired 
tolerances.  The  mole  fractions  of  the  two  phases  are  then  given  by 

*  exp  (ziQ)  (A.  14) 

y.  *  exp  (w±o)  (A. 15) 

The  computer  program  is  listed  as  shown  on  the  following 

pages . 


IMPLICIT  REAL*8( A-Hf  0-Z ) 

PREDICTION  OF  TERNARY  MIXTURE  IMMI SCIBI LI TY 
VAN  LAAR  MODEL  UTILIZED 
STOICHIOMETRIC  REQUIREMENT  INCLUDED 

DIMENSION  Z(  12)  ,BU2),C(  12)  »D112),E(  12)  ,F(  12)  ,G(  12) ,DHF< 12)  ,H(12) 
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APPENDIX  B 


THERMO CHEMICAL  PROPERTY  DATA 

B.l  Thermodynamic  Properties  of  Li(s,l,g),  Li2(g),  LiF(s,l,g), 

L1JfJ(I)7  Li3F3(g),  and  SF6<g) 

The  thermodynamic  properties  of  Li(s,l,g),  Li2(g),  LiF(s,l,g), 
Li2F2(g),  Li^F^(g)  and  SF^(g)  were  taken  from  the  JANAF  Thermo chemical 
Tables  [10].  The  symbols  s,  1,  and  g  signify  solid,  liquid,  and 
gas  phases,  respectively.  The  data  were  inserted  into  the  computer 
program  by  utilizing  polynominal  curve  fits  available  from  JANAF. 

B.2  Thermodynamic  Properties  of  Li2S(s,l,g) 

The  information  available  in  the  literature  on  the  properties 
of  Li2S  is  very  limited.  The  heat  of  formation  of  Li2S(s)  at  298°K 
was  calculated  by  Morris  [13]  to  be  -492.4+21  kJ/mole,  and  determined 
experimentally  by  Juza  and  Uphoff  [14]  to  be  -445.6  kJ/mole.  The 
melting  point  of  Li2S  was  experimentally  determined  to  be  1645+^°K 
by  Cunningham,  Johnson,  and  Cairns  [15],  The  boiling  point  of  Li2S 
was  determined  to  be  2073°K  by  Mott  [16]  in  an  experimental  study 
of  the  volatilities  of  various  refractory  materials.  Since  the 
above  information  was  incomplete,  it  was  necessary  to  estimate  the 
remaining  thermodynamic  property  data. 

B.2.1  Thermodynamic  Properties  of  Li^SCs).  The  specific  heats 
of  Li2S(s)  were  estimated  by  using  Knopp's  rule  [38]  with  data  taken 
from  Reference  [10]  for  Na2S(s),  Li20(s)  and  Na20(s).  The 
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estimated  heat  capacity  data  were  then  fitted  to  a  polynominal 
equation  using  least  squares  methods.  The  estimated  heat  capacities 
were  found  to  be  approximately  7  to  20%  lower  than  the  corresponding 
values  for  Na^SCs)  depending  on  the  temperature.  Estimates  of  the 
entropy  of  formation  of  L^SCs)  were  given  by  Steunenberg  [47]  at 
temperatures  of  623,  673,  and  723°K.  The  estimated  entropies  of 
formation  were  used  to  determine  the  entropies  of  Li^SCs)  at  the 
given  temperatures,  which  were  in  turn  used  in  conjunction  with  the 
estimated  heat  capacity  data  to  determine  the  entropy  at  298°K. 

The  enthalpies  of  Li^SCs)  were  determined  by  integration  of  the  heat 
capacity  equation  and  using  the  value  of  Morris  [13]  for  the  heat 
of  formation  at  298°K. 

B.2.2  Thermodynamic  Properties  of  Li2S(l).  The  heat 
capacities  of  Li2S(l)  were  estimated  by  using  Knoop’s  rule  with  data 
taken  from  Reference  [10]  for  Na2S(l),  Li20(l)  and  Na20(l).  The 
estimated  data  were  found  to  be  4.2  J/mole-°K  lower  than  the 
corresponding  values  for  Na2S(l)  and  the  data  were  represented  by  the 
JANAF  polynominal  equation  for  the  Na2S(l)  data  after  the  proper 
adjustment.  The  heat  of  fusion  of  lithium  sulfide  was  estimated  by 
compiling  data  from  References  [10]  and  [48]  through  [51]  on  the 
melting  points  and  heats  of  fusion  of  alkali  metal-halogen  salts, 
lithium  and  sodium  oxides,  alkali  metal  sulfates,  and  sodium  sulfide. 

i 

! 

The  entropies  of  fusion,  ASf,  were  obtained  by  noting: 


(B.l) 
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The  entropies  of  fusion  of  the  various  compounds  were  plotted 
against  the  molecular  weight  of  the  participating  alkali  metal. 

Definite  trends  were  noticed  in  the  curves,  especially  between  the 
lithium  and  sodium  compounds.  By  continuing  the  trends,  the  entropy 
of  fusion  of  Li2S  was  estimated  to  be  3. 8+1. 2  J/mole-°K.  By 
employing  Equation  (B.l)  and  the  melting  point  of  Cunningham,  et  al. , 
[15],  the  heat  of  fusion  was  estimated  to  be  6.19+2.1  kJ/mole.  The 
estimated  value  was  taken  to  be  7.24  kJ/mole  in  subsequent  calculations 
as  this  value  had  been  used  in  earlier  work.  A  value  for  the  heat  of 
fusion  of  49.4  kJ/mole  was  also  considered  since  this  value  resulted 
in  the  best  representation  of  the  lithium  fluoride-lithium  sulfide 
binary  solubility  data.  The  heat  of  formation  of  a  solid  and  liquid 
can  be  related 
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where  AHf  is  the  heat  of  fusion  at  the  melting  point  of  the  pure 
component  and  AHfQrm  is  the  heat  of  formation.  The  variation  of  the 
heat  of  formation  with  temperature  can  be  expressed  as  follows: 
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Combination  of  Equations  (B.2)  and  (B.3)  results  in  an 
expression  relating  the  heat  of  formation  of  Li^Sd)  at  298°K  to  the 
heat  of  formation  of  L^SCs)  at  298°K,  the  heat  of  fusion,  and 
relative  enthalpies  of  the  compounds  as  follows: 
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(B.4) 


Utilizing  Equation  (B.4),  the  heat  of  formation  of  L^SCl)  was 
estimated  to  be  -493.3  kJ/mole  for  the  low  heat  of  fusion  case  and 
-451.0  kJ/mole  for  the  high  heat  of  fusion  case.  The  entropy  of 
L^SCl)  at  298°K  was  equated  to  the  entropy  of  Li^SCs)  at  298°K 
plus  the  difference  between  the  corresponding  entropies  for  Na^SQ) 
and  Na^SCs).  The  enthalpies  of  Li.Sd)  were  determined  by  integration 
of  the  heat  capacity  equation  and  using  the  above  values  for  the 
heats  of  formation  at  298°K. 


B.2.3  Thermodynamic  Properties  of  L^SCg).  The 
thermodynamic  properties  of  Li^SCg)  were  estimated  by  application  of 
the  principles  of  statistical  thermodynamics  after  the  initial 
estimation  of  the  required  spectroscopic  data. 

The  vibrational  frequencies  of  the  molecule  were  estimated  by 
a  comparison  of  the  trends  in  similar  diatomic  and  triatomic  molecules. 
Data  taken  from  the  JANAF  Thermo chemical  Tables  [10]  and  used  for 
this  purpose  are  listed  in  Table  6,  The  vibrational  frequencies 
were  estimated  simply  by  continuing  the  trends  shown  in  molecules 
having  hydrogen  atoms  substituted  for  lithium  atoms.  The  first 
step  was  to  estimate  the  simple  harmonic  oscillator  frequency  of  LiS 
from  trends  shown  in  HN,  HO,  HS,  LiN  and  LiO.  It  was  then  observed 
that  the  frequency  of  the  triatomic  v  mode  was  approximately  3% 
larger  than  the  corresponding  diatomic  harmonic  oscillator  frequency 
for  the  ^N,  ^0,  H^S  and  Li^O  molecules.  The  frequency  of  Li^S 
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was  estimated  accordingly.  The  and  frequencies  of  the  bent 
molecules  were  observed  to  be  approximately  in  similar  proportion  to 
the  corresponding  frequencies.  As  a  result  the  and 
frequencies  of  L^S  were  taken  to  be  in  the  same  proportion  to  the 
frequency  as  in  the  Li^O  case. 

This  method  of  estimation  of  the  vibrational  frequencies  was 
expected  to  give  only  relatively  crude  estimates.  The  effects  of 
changing  these  values  on  the  thermodynamic  properties  was  examined. 

The  Li-S  bond  distance  was  estimated  from  the  Li-0  bond 
distance  in  Li^O  by  assuming  the  same  fraction  of  the  sum  of  the 
ionic  bond  radii  of  the  participating  atoms.  The  atomic  radii  were 
taken  from  Lange  [48],  For  L^O,  the  Li-0  bond  distance  was  found 
to  equal  the  sum  of  the  atomic  radii  multiplied  by  a  factor  of  0.779. 
Hence  the  Li-S  bond  distance  in  Li2S  was  estimated  to  be  1.93  A.  The 
rotational  constants  were  then  computed  using  this  estimated  value 
of  the  bond  distance.  The  accuracy  of  this  method  was  examined  by 
applying  the  same  method  to  the  ^0 ,  and  HjS  molecules .  The 
accuracy  depends  on  which  molecule  was  picked  to  determine  the  ratio, 
but  the  maximum  error  that  resulted  was  only  2.8%.  The  effects  of 
varying  the  bond  distance  was  also  examined. 

The  statistical  thermodynamic  equations  utilized  to  compute 
the  thermodynamic  properties  are  summarized  below.  The  equations  are 
developed  in  any  statistical  thermodynamics  text  such  as  the  text  by 
Davidson  [52],  All  units  are  in  J-mole  ^  or  J-mole  '*'-0K 

The  electronic  ground  state  of  L^S  was  assumed  to  be  a  singlet 
with  no  low-lying  electronic  states.  This  is  the  case  for  L^O.  Hence 
there  were  no  electronic  contributions  to  the  thermodynamic  properties. 
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The  translational  contributions  to  the  thermodynamic  functions 

were  computed  by  the  following  for  both  the  linear  and  nonlinear 

_2 

olecules.  The  reference  pressure  was  taken  to  be  101325  n-m  , 

-3 

and  Mq  equals  1  x  10  kg. 
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The  vibrational  contributions  are  given  by  the  following: 
v.(cm-1) 

where  u^^  -  1.4388  ■  — 
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The  bending  vibration  is  nondegenerate,  ■  1,  for  the  bent 
molecule  and  doubly  degenerate,  ■  2,  for  the  linear  molecule. 


(B. 11) 
(B. 12) 
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The  symmetry  number  (a)  for  both  the  linear  and  bent  molecules 

* 

is  equal  to  two  due  to  the  symmetric  structure  of  the  molecule. 

The  total  thermodynamic  quantities  are  given  by  summation  of 
the  electronic,  translational,  rotational  and  vibrational 
contributions.  The  physical  constants  used  in  the  calculation  were 
taken  from  Davidson  [52]. 

Properties  were  computed  for  both  the  linear  molecule  and 
bent  molecule  utilizing  the  estimated  data.  The  Li-S-Li  bond 
angle  was  taken  to  be  100  degrees  in  the  bent  case.  To  examine 
changes  brought  about  by  variances  in  the  assumed  data,  four  additional 
cases  were  examined  for  each  configuration  as  follows :  (a)  vibrational 

frequencies  doubled,  bond  distance  fixed;  (b)  vibrational  frequencies 
halved,  bond  distance  fixed;  (c)  base  vibrational  frequencies,  bond 
distance  reduced  by  10%;  and  (d)  base  vibrational  frequencies,  bond 
distance  increased  by  10%. 

Table  7  lists  the  percent  deviation  of  the  indicated 
thermodynamic  quantities  from  the  values  given  by  the  base  case  at 
a  temperature  of  1200°K„  The  base  case  was  the  linear  molecule  in 
conjunction  with  previously  estimated  spectroscopic  data.  It  is 
believed  that  the  estimated  spectroscopic  data  is  much  more  precise 
than  the  variances  examined  and  therefore  it  is  believed  that  the 
estimated  thermodynamic  data  is  accurate  to  within  +5%. 

The  heat  of  formation  of  L^SCg)  was  required  to  determine  the 
values  of  the  absolute  enthalpies.  The  heat  of  formation  was 
determined  by  choosing  a  value  which  when  used  in  conjunction  with 
the  data  previously  estimated,  computed  the  boiling  point  given  by 
Mott  [16].  The  values  determined  were  -220.2  kJ/mole  for  the  low 
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Table  7 


Effects  of  Various  Assumptions  on  the  Thermodynamic  Properties 


Case 

0 

c 

Linear** 

Bent 

s° 

Linear 

Bent 

-(f°-ho 

Linear 

°/T) 

Bent 

h°-h 

Linear 

o 

o„ 

Bent 

Base  Case 

0.0 

-6.7 

0.0 

+2.4 

0.0 

+4.2 

0.0 

-6.6 

Vibrational 

Frequency 

Doubled 

-3.6 

-10.2 

-6.7 

-2.6 

-6.2 

-0.1 

-8.9 

-14.6 

Vibrational 

Frequency 

Halved 

-1.0 

-5.8 

+6.9 

+7.6 

+7.3 

+9.4 

+5.5 

-1.5 

Bond  Distance 
Decreased  102 

0.0 

-6.7 

-0.5 

+1.6 

-0.6 

+3.2 

0.0 

-6.6 

Bond  Distance 
Increased  102 

0.0 

-6.7 

+0.5 

+3.1 

+0.6 

+5.1 

0.0 

-6.6 
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heat  of  fusion  case,  and  -177.9  kJ/mole  for  the  high  heat  of  fusion 
case. 

A  check  on  the  above  estimation  was  made  by  computation  of  the 
implied  heat  of  vaporization.  At  the  boiling  point  the  following 
holds: 


41ib  '  “fora,  <«>  '  “fora,  (1>  <»•»> 

b  b 

where  AH^  is  the  heat  of  vaporization  and  T^  is  the  boiling  temperature. 
Substitution  of  Equation  (B.3)  yields: 


The  value  for  the  heat  of  vaporization  was  determined  to  be 
230  kJ/mole  for  both  the  high  and  low  heats  of  fusion  cases.  The 
value  agrees  well  with  the  value  of  182  kJ/mole  determined  by 
application  of  Trouton’s  rule  [38]  which  states  that  the  entropy 
increase  upon  boiling  is  87.9  J/mole-°K. 


B.3  Densities  of  Li(l),  LiF(l)  and  Li2S(l) 

The  density  of  pure  lithium  was  taken  from  Davison  [12],  and 
the  density  of  pure  lithium  fluoride  was  taken  from  Janz  [11].  The 
density  of  pure  liquid  lithium  sulfide  is  not  reported  in  the 
literature.  The  density  of  solid  lithium  sulfide  at  298°K  was  given 

3 

by  Mourlot  [17]  to  be  in  the  range  1630  to  1700  kg/m  .  Zintl,  et  al., 
[18]  examined  the  structure  of  the  solid  and  determined  the  density 
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-i 

to  be  1640  kg/m  at  298°K.  Beck  [19]  presented  a  correlation 

between  the  volume  contraction  of  the  constituent  elements  and  the 

energy  of  formation  upon  chemical  reaction.  The  correlation  yielded 

3 

a  value  for  the  density  of  solid  lithium  sulfide  of  1640  kg/m  . 

3 

The  density  at  298°K  was  selected  to  be  1660  kg/m  ,  and  the  density 

variation  with  temperature  was  estimated  to  have  a  slope  equal  to  the 

average  slope  of  lithium  and  lithium  fluoride.  Density  changes  upon 

melting  were  neglected  since  such  changes  are  small.  The  equations 

3 

representing  the  variation  of  the  density  in  kg/m  of  the  various 
components  with  temperature  are 


l  .  =  562.  -  0 . IT 

Li 

(B. 25) 

>T  =  2358,1  -  0.4902T 

Lir 

(B . 26) 

>TJ  =  1748.  -  0.294T 

(B. 27) 

Li2S 

APPENDIX  C 


THEORETICAL  DEVELOPMENT  OF  EXPRESSIONS  USED  TO  DETERMINE 
THERMODYNAMIC  MODEL  PARAMETERS 


C. 1  Expressions  Relating  to  Solid-Liquid  Equilibria 

The  development  presented  below  is  based  on  a  similar  analysis 

given  in  Reference  [34].  At  the  freezing  temperature,  Tf  ,  of 

i 

component  i,  the  free  energies  of  the  pure  solid  and  pure  liquid 
are  related  as  follows: 


f±(s) 


(C.l) 


If  it  is  assumed  that  solid  solutions  are  not  formed,  then  the  free 
energies  of  the  solid  and  liquid  phases  in  equilibrium  at  a  temperature 
T  and  a  mole  fraction  are  related  by 


ffCs)  T  *  f^solu,  x^  T 


(C.  2) 


Equations  (C.l)  and  (C.2)  are  combined  and  the  effects  of  composition 
and  temperature  changes  are  separated  by  writing 


(f1(solu,  xi)~fiC 


T  -  (fi(s)-fi(solu,  xi))|T 
f  i 


-  (fi(s)-fi(solu,  xi>) 


(C.  3) 


Equation  (2.5)  is  evaluated  at  temperature  T,  ,  and  is  substituted  into 

i 

Equation  (C.3).  The  resulting  equation  is  rearranged  and  expressed 
as  the  following  Integral  relationship: 
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R  In  a. 


r  f^s)-! i(solu,  x^) 

“  /  d  - T - 

C1  JT *  L  J 


(C.4) 


It  Is  noted  that  at  constant  pressure 


a||i)  -  ^  ai 


(C.5) 


and  Equation  (C.5)  can  be  written  as 

T 


R  In  a. 


■/ 


(solu,  xi)-hi(s) 


dT 


i  T, 


(C.6) 


The  relative  molal  heat  content  is  defined  as 


(C.  7) 


(C.8) 


ki  S  h^Csolu,  x^-hj0 

with  the  molal  enthalpy  h^solu,  x^  defined  as 


i-  ,  .  3H 

h^SOlU,  Xt)  - 

uj(j  *  i),  P,  T 

The  heat  of  fusion  is  approximated  by 

T 

AH  -  h  °  -  h. (s)  -  AH,  +  Ac  (T  -  T,  ) 
fi  1  1  fi  pi  fi 

Tf 

where  AH  is  the  standard  heat  of  fusion  of  component  i  at  its 
*i 

melting  temperature,  and  Ac^  is  the  difference  between  the  partial 
molal  heat  capacity  of  component  i  in  solution  and  the  molal  heat 
capacity  of  the  solid.  Substitution  of  Equations  (C.7)  and  (C.9) 
into  Equation  (C.6)  yields 


(C.9) 
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R  In  a. 


(C.10) 


where 


e  -  t  -  t 

i 


(C.ll) 


The  right-hand-side  of  Equation  (C.10)  can  be  expressed  as  the  sum  of 
two  parts:  the  term  within  the  integral  is  the  contribution  from  the 
enthalpy  change  of  nonideal  mixing,  and  the  remaining  terms  are  the 
contributions  from  the  heat  of  fusion  of  component  i. 

The  van  Laar  thermodynamic  model,  Equation  (2.10),  is  used  to 
evaluate  both  the  relative  molal  heat  content  and  the  activity  of 
component  i.  Equations  (2.20)  are  evaluated  at  temperature  T 

ri 

to  yield  the  left-hand-side  of  Equation  (C.10).  The  relative  molal 
heat  content  expressions  are  developed  from  Equation  (2.30).  It  is 
assumed  that  the  empirical  parameters  representing  the  binary  system, 
and  B^,  vary  linearly  with  temperature  as  follows: 

A±j  -  A  +  BT 

Bi  -  a  +  bT  (C.12) 


Employing  Equations  (C.12),  the  expressions  for  the  relative  heat 
contents  of  components  1  and  2  in  Binary  1-2  are 


(C. 13) 
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Substitution  of  the  expression  for  the  activity  of  component  1  given 
by  Equation  (2.20)  and  Equation  (C.13)  into  Equation  (C.10)  yields 
an  equation  representing  the  thermodynamic  equilibrium  between  solid 
component  1,  and  component  1  in  solution.  Substitution  of  the 
expression  for  the  activity  of  component  2  given  by  Equation  (2.20) 
and  Equation  (C.14)  into  Equation  (C.10)  yields  a  similar  expression 
for  component  2. 

Equations  (C.10)  can  be  solved  for  the  liquid  composition  as 
a  function  of  temperature  if  the  van  Laar  empirical  parameters  are 
known:  conversely,  the  equations  can  be  used  in  conjunction  with 
known  phase  equilibria  data  to  compute  the  empirical  parameters. 

The  most  common  phase  equilibria  data  employed  with  the  equations 
are  determined  from  the  liquidus  line  or  freezing  point  depression 
line  obtained  from  the  lowering  of  the  freezing  point  of  the  pure 
solvent  on  the  addition  of  solute. 


C. 2  Expressions  Relating  to  Two  Equilibrated  Liquid  Phases 

As  discussed  in  Section  2.3,  the  requirement  for  two  phases 
to  exist  in  equilibrium  reduces  to  the  activities  of  each  component 
being  identical  in  each  of  the  two  phases.  If  the  two  phases  are 
indicated  by  mole  fractions  x^  and  y^,  the  requirement  can  be  stated 
as 


ai(xi) 


ai(yi) 


(C. 15) 
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The  van  Laar  model  for  the  activity  coefficients,  Equations  (2.20), 


is  employed  to  yield 
A 
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(C. 16) 


Equations  (C.16)  can  be  solved  for  the  two  unknown  empirical 
parameters,  and  A^2,  as  a  function  of  temperature  and  the 
equilibrium  compositions  as  follows: 
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(C. 18) 


Equations  (C.17)  and  (C.18)  are  solved  for  the  empirical 
parameters  at  a  given  temperature  from  knowledge  of  the  compositions 
of  the  two  equilibrated  liquid  phases  at  that  temperature. 
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BATH  SAMPLE  ANALYSIS  PROCEDURE 

D. 1  Lithium  Fluoride  and  Lithium  Sulfide  Determination 

The  number  of  moles  of  lithium  fluoride  and  lithium  sulfide 
in  both  the  fuel-rich  and  product-rich  bath  samples  were  determined 
using  the  following  wet-test  chemical  analysis  procedure. 

D.1.1  Sulfide  Separation.  The  lithium  sulfide  was  separated 
from  the  samples  using  a  distillation  procedure  described  by  Luke  [53]. 
A  schematic  of  the  apparatus  is  shown  in  Figure  25,  A  500  ml 
Erlenmeyer  flask  was  used  for  both  the  sample  flask  and  the  receiving 
flask. 

The  sample  which  had  been  prepared  for  analysis  was  introduced 
to  the  sample  flask  under  an  argon  atmosphere  in  a  dry  box.  The 
sample  flask  was  then  introduced  to  the  distillation  apparatus. 

A  sufficient  quantity  of  distilled  water  was  added  through  the 
separatory  funnel  to  ensure  complete  solution  of  the  lithium  fluoride. 
The  solubility  of  lithium  fluoride  in  water  is  approximately  0,13 

3 

gm  per  100  cm  at  298CK.  For  fuel-rich  samples  of  approximately  1.0 

3 

gm  and  product-rich  sample  of  approximately  0.10  gm,  300  and  350  cm 

of  water  were  added,  respectively.  The  sample  solution  was  acidified 
3 

by  adding  130  cm  of  one  molar  hydrochloric  acid  to  the  fuel-rich 

3 

samples  and  50  cm  of  half  molar  hydrochloric  acid  to  the  product- 


rich  samples.  An  argon  purge  transferred  the  evolved  hydrogen 
sulfide  gas  to  the  receiving  flask  positioned  in  an  ice  bath  and 
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3  3 

containing  350  cm  of  water  and  40  cm  of  ammonical  zinc  sulfate 
solution.  The  sample  solution  was  heated  for  a  total  of  36  minutes 
of  which  approximately  12  minutes  were  at  a  slow  boil  condition. 

The  sample  flask  contained  F  and  Li+  ions  from  the  three  sources  of 
free  lithium,  lithium  fluoride  and  lithium  sulfide. 

D.1.2  Lithium  Sulfide  Determination.  Reagents: 

1.  Ammonical  Zinc  Sulfate  Solution  -  dissolve  50.0  gm 

ZnS04. 7H  0  in  250  ml  NH  OH  (28%  NHj. 

2  4  3 

2.  Potassium  Iodide  -  KI  -  0.50  M. 

3.  Potassium  Iodate  -  KIO^  -  0.05  M. 

4.  Soluble  Starch  -  1%  in  H^O. 

The  procedure  described  by  Luke  [53]  was  utilized.  The 
solution  in  the  receiving  flask  was  transferred  to  a  larger  flask  with 
the  aid  of  distilled  water.  Sufficient  Kl  and  HC1  was  added  to 
provide  an  excess  of  H  and  I  in  the  solution.  Starch  indicator 
was  added,  and  the  solution  was  titrated  with  KIO^  until  the 
appearance  of  a  permanent  blue  color  indicating  the  presence  of  I^. 

The  following  equations  described  the  process: 

10 3“  +  5I~  +  6H+  -*•  3I2  +  3H20  (D.l) 

H2S  +  I2  2H1  +  s  (D.2) 

The  number  of  moles  of  H2S  present  is  three  times  the  number 
of  moles  of  10^  added  during  the  titration.  The  H2S  was  produced 
by  the  following  reaction: 


Li2S  +  H20  -*•  2L10H  +  HgS 


(D.3) 
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Therefore  the  number  of  moles  of  Ll^S  present  also  equals 
three  times  the  number  of  moles  of  10  ^  added. 

D.1.3  Lithium  Fluoride  Determination.  Reagents: 

1.  Thorium  Nitrate  -  TH(No»)  -  0.025  M. 

J  4 

The  procedure  described  by  Willard  and  Winter  [54]  was 

utilized.  The  solution  remaining  in  the  sample  flask  was  transferred 

3 

to  a  graduated  cylinder  and  the  volume  determined.  A  20  cm  aliquot 
was  measured  with  a  hypodermic  syringe  and  transferred  to  a  graduated 

3 

cylinder.  Twenty  cm  of  ethyl  alcohol  were  added  to  the  cylinder 

along  with  three  drops  of  alizarin  Red  S  indicator  (1%).  The 

solution  was  titrated  with  thorium  nitrate  until  the  appearance  of  a 

faint  pink  color.  The  thorium  nitrate  precipitated  thorium  fluoride 

which  is  insoluble  in  ethyl  alcohol,  and  the  indicator  detected  the 

point  where  all  the  fluoride  was  precipitated.  The  titration  was 

3 

repeated  for  five  twenty  cm  samples.  The  following  equation 
describes  the  process: 

Th+4  +  4  F~  ThF4  (D.4) 

The  number  of  moles  of  F  present  is  four  times  the  number  of  moles 

of  Th (NO  )  added  multiplied  by  the  ratio  of  the  total  solution 
4 

volume  to  the  sample  volume. 

D.2  Lithium  Determination 

Different  methods  were  used  to  determine  the  lithium  content 
of  the  product-rich  and  fuel-rich  samples . 
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D.2.1  Product-Rich  Samples.  The  method  of  determining  the 
free  lithium  contained  in  the  product-rich  sample  was  based  on 
the  measurement  of  the  released  by  the  following  reaction. 

Li  +  H20  +  LiOH  +  1/2  H2  (D.5) 

The  combustor  sample  (approximately  10  grams)  was  accurately 

3 

weighed  and  placed  in  a  glass  chamber  with  a  volume  of  9957  cm  .  The 

chamber  was  evacuated  and  backfilled  with  argon  gas  to  atmospheric 

pressure.  The  chamber  pressure  was  measured  with  a  Meriam  Model 

31EC10,  60  inch,  mercury  filled  absolute  manometer.  A  known  quantity 

3 

of  distilled  water  (approximately  450  cm  )  was  added  to  the  chamber. 
After  twenty  minutes,  the  system  pressure  and  temperature  were 
recorded.  The  total  number  of  moles  of  argon,  water  vapor,  hydrogen 
sulfide  and  hydrogen  present  in  the  sample  chamber  were  therefore 
known  by  the  ideal  gas  law. 

The  hydrogen  measurement  system  consisted  of  a  Fisher-Hamilton 
Model  29  gas  partitioner  equipped  with  a  80/100  mesh  5-A  molecular 
sieve.  The  gas  sample  was  introduced  to  the  partitioner  through  a 

3 

0.5  cm  sample  loop  which  is  integral  with  the  device.  The  output 
from  the  thermal  conductivity  detector  was  recorded  on  a  Leeds  and 
Northrup  Model  W  Speedomax  recorder.  The  detection  system  was 
calibrated  with  prepared  samples  of  H2  in  argon.  Therefore,  from 
the  recorder  output  the  ratio  of  the  number  of  moles  of  H2  to  the 
number  of  moles  of  argon  in  the  sample  was  determined. 

The  gas  partitioner  was  connected  to  the  sample  chamber. 

The  gas  sample  circuit  was  evacuated  and  then  opened  to  the  sample 

3 

chamber.  A  gas  sample  of  approximately  80  cm  was  extracted  from 
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the  sample  chamber.  The  gas  sample  pressure  was  reduced  to  atmospheric 
pressure  and  sent  through  the  gas  partitioner  column  for  detection. 

The  output  gave  the  ratio  of  the  number  of  moles  of  to  the  number 
of  moles  of  argon,  water  vapor,  and  hydrogen  sulfide.  The  gas  sample 
pressure  and  temperature  had  to  be  corrected  to  that  of  the  calibration 
standard.  The  information  from  the  gas  partitioner  combined  with 
the  knowledge  of  the  total  number  of  moles  of  gas  present  in  the  sample 
chamber  yielded  the  total  number  of  moles  of  evolved.  Employing 
Equation  (D.5)  yielded  the  total  number  of  moles  of  free  lithium  in 
the  sample. 


D.2.2  Fuel-Rich  Samples.  The  free  lithium  content  of  the  fuel- 
rich  samples  was  obtained  by  subtracting  the  sum  of  the  lithium 
fluoride  and  lithium  sulfide  weights  from  the  total  sample  weight. 

The  difference  procedure  was  accurate  for  the  samples  analyzed 
because  the  samples  were  greater  than  95%  lithium  on  a  mole  fraction 
basis .  A  procedure  was  examined  that  permitted  the  direct  measurement 
of  the  total  Li+  content  of  the  sample.  The  procedure  can  be  applied 
in  situations  where  the  difference  method  cannot  be  used  accurately. 

The  procedure  is  described  by  Rogers  and  Calley  [55], 

Reagents: 

1.  Periodate  Reagent  -  24.0  gm  KOH,  10.0  gm  KIO^,  100  ml  I^O. 

2.  Sodium  Thiosulfate  -  Na2®2°3  ~ 

3.  Sulfuric  Acid  -  H^O^  -  1.0  N. 

4.  Potassium  Hydroxide  -  KOH  -  4.0  N. 

3 

A  5.0  cm  aliquot  of  the  sample  solution  was  measured  with  a 

3 

hypodermic  syringe  and  transferred  to  a  50  cm  beaker  maintained  at 
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60°-70°C  by  a  heated  water  bath.  The  periodate  reagent  was  added 
in  a  dropwise  manner  from  a  buret,  and  the  lithium  was  precipitated 
as  a  complex  periodate  mixture  (y  Li^IO^).  The  precipitate  was 
filtered  and  washed  with  KOH,  and  dissolved  in  H.SO. .  Sufficient  KI 
was  added  to  ensure  complete  conversion  of  the  lithium  periodate  to 
I 2*  Starch  indicator  was  added  to  indicate  the  presence  of 
The  solution  was  titrated  with  Na^S^O^  until  the  disappearance  of  the 
blue  color.  The  following  equations  describe  the  process: 

10  ~5  +  11  +  12H+  -*•  41.  +  6H-0  (D.  6) 

O  l  L 

2Na2S2°3  +  I2  +  2NaI  +  Na2S4°6  (D,7) 

The  number  of  moles  of  Li  present  is  equal  to  5/8  of  the  number  of 
moles  of  Na2S2°3  titrated  to  the  solution  multiplied  by  the  ratio  of 
the  total  solution  volume  to  the  sample  volume.  The  process  was 
calibrated  with  known  lithium  samples  for  a  range  of  lithium 
concentrations,  since  the  ratio  can  vary  slightly  with  concentrations. 
The  above  analysis  was  repeated  for  five  samples. 
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APPENDIX  E 

OXIDIZER  UTILIZATION  TESTS 

E. I  Apparatus  and  Procedure 

The  oxidizer  utilization  efficiency  tests  were  conducted 
on  two  radiatively  cooled  pot-type  combustors  of  different  size. 

The  oxidizer  was  injected  into  the  molten  bath  through  a  coaxial 
injector  mounted  at  depths  of  5.0  to  7.6  cm  below  the  fuel  surface. 

The  utilization  efficiency  was  determined  by  measuring  the  oxidizer 
concentration  in  the  argon  gas  exhausted  from  the  combustor  at  a  known 
rate.  The  oxidizer  flow  rate  from  the  combustor  was  then  compared  to 
the  oxidizer  flow  rate  injected  into  the  bath. 

3 

Samples  of  the  exhaust  flow  were  collected  in  125  cm 
evacuated  sampling  bulbs  manufactured  by  the  Arthur  H.  Thomas  Company 
(P/N  5585-K15) .  The  lithium  vapor  was  removed  from  the  sample  prior 
to  entry  into  the  sampling  bulb  to  protect  the  gas  analysis  system 
from  contamination  by  passing  the  sample  through  a  one-pass  felt 
metal  screen.  Calibration  tests  revealed  no  measurable  influence  of 
the  screen  on  the  oxidizer  concentration  in  the  sample. 

A  Perkin-Elmer  model  880  gas  chromatograph  was  used  for  the 
analysis  of  the  argon-oxidizer  gas  mixture.  The  separating  column 
was  stainless  steel,  6.35  mm  O.D.  by  122  cm  in  length,  and  was 
packed  with  80-100  mesh  silica  gel.  The  column  was  placed  in  an 
ice  bath  and  operated  at  0°C.  The  system  was  calibrated  with  prepared 
samples  of  SF^-argon  mixtures.  The  prepared  samples  had  concentration 
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of  SF^  of  938,  1930,  5680  and  7050  parts  per  million.  Preparation 
tolerances  were  estimated  to  be  ±10%  for  the  low  concentrations  and 
±3%  for  the  high  concentrations.  A  commercially  prepared  sample  of 
concentration  106  parts  per  million  with  a  tolerance  of  ±2%  was 
purchased  from  Matheson  Gas  Products.  The  calibration  procedure 
indicated  that  the  SF^  concentration  in  the  argon  was  a  linear 
function  of  the  SF^  peak  area  in  the  concentration  range  of  interest. 
Helium  was  used  as  the  carrier  gas .  SF,  concentrations  as  low  as  10 

D 

parts  per  million  could  be  detected  with  the  apparatus. 

E.2  Results 

Samples  were  taken  approximately  every  10  to  20  minutes  during 
the  1.5  to  2  hour  period  following  ignition.  No  oxidizer  could  be 
detected  in  the  exhaust  flow  after  the  second  sample  taken  at  the  20 
minute  mark.  Bath  temperatures  at  this  time  were  in  the  range  of  944 
to  1061°K  with  ignition  occurring  at  a  bath  temperature  of  572°K. 

Small  quantities  of  oxidizer  were  detected  in  samples  taken  during 
the  first  15  minutes  following  ignition,  but  the  oxidizer  utilization 
efficiency  exceeded  99.6%  for  all  measurements.  It  was  concluded 
from  the  tests  that  the  oxidizer  utilization  efficiency  is  100%  during 
normal  combustor  operation  with  Injector  submergence  depths  exceeding 
5  cm. 
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APPARATUS  AND  PROCEDURE  FOR  INVESTIGATING  THE 
WET-TANK  PRODUCT  COLLECTION  METHOD 


F.l  Apparatus 

A  schematic  of  the  apparatus  used  to  investigate  the 

feasibility  of  collecting  and  storing  the  reaction  products  in  the 

fuel  tank  is  shown  in  Figure  26.  The  product  liquid  was  maintained 

at  a  temperature  of  1227 °K  in  the  upper  vessel,  and  the  lithium  fuel 

was  maintained  at  a  temperature  of  505° K  in  the  lower  vessel. 

The  temperature  of  the  vessels  were  maintained  by  heaters  independently 

controlled  with  variacs.  The  product  liquid  was  forced  from  the 

3 

upper  vessel  at  a  flow  rate  of  35-55  cm  /min.  into  the  lower  vessel 
by  an  inflow  of  argon  gas  into  the  upper  vessel.  Temperatures  were 
measured  with  chromel-alumel  thermocouples  and  recorded  on  a  Leeds 
and  Northrup  Speedomax-H  null-balance  recorder. 

F.2  Procedure 

The  apparatus  was  stabilized  at  the  desired  temperature  levels 
with  an  argon  purge  supplied  to  the  upper  vessel.  The  purge  flow 
was  exhausted  through  vents  in  the  upper  and  lower  vessels.  Product 
liquid  outflow  was  initiated  by  isolating  the  ullage  region  of  the 
upper  vessel  using  a  series  of  valves  and  pressurizing  the  upper 
vessel  with  the  argon  inflow.  The  flow  rate  was  metered  by  metering 
the  argon  Inflow  and  correcting  the  rate  for  thermal  expansion  in  the 
vessel.  An  additional  estimate  of  the  flow  rate  was  obtained  by 
estimating  the  outflow  period  from  the  temperature-time  curve  of  the 
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Figure  26  Wet-Tank  Product  Collection  Apparatus 


fuel  bath  and  knowing  the  mass  and  density  of  the  material 
exhausted.  For  tests  where  pure  lithium  fluoride  was  used  to 
simulate  the  reaction  products,  the  form  of  the  exhausted  material 
was  observed  by  removing  the  lithium  from  the  lower  vessel.  The 
lithium  was  removed  by  reacting  it  with  water  after  the  system  had 
cooled.  For  tests  where  actual  reaction  products  were  used,  the 
lithium  could  not  be  removed  with  water  because  the  lithium  sulfide 
component  is  soluble  in  water  and  the  form  of  the  exhausted 
products  would  be  altered.  The  exhausted  product  was  retrieved  from 
the  molten  fuel  bath  with  a  stainless  steel  mesh  basket  placed  in  the 
bottom  of  the  lower  vessel. 

F.3  Results 

The  results  indicated  that  product  liquid  exhausted  into 
a  fuel  tank  at  rates  typical  of  the  steady  combustor  examined  in  this 
study  would  form  small  spherical  balls  approximately  0.86  cm  in 
diameter  with  a  hollow  core.  Identical  balls  were  found  when  either 
pure  lithium  fluoride  or  actual  reaction  products  were  used  in  the 
tests.  The  packing  efficiency  of  the  balls  was  found  to  be 
approximately  50%  on  a  volume  basis. 
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